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Abstract 

Applications of QCD sum rules to exclusive decays of heavy mesons are reviewed. 
In detail, we discuss the calculations of the B and D decay constants, the heavy- 
Q\ ' to-light form factors f + (p 2 ) and f°{p 2 ), and the strong couplings gB*Bir and qd*d-k- 

Predictions are presented for the semileptonic weak decays B — > idv\ and B — > plvi, 
and used to extract V u b from the measured decay widths. Results are also shown on 
| D — > ixlvi and D* — ► L>7r. Furthermore, we reconsider the factorization hypothesis in 

nonleptonic two-body decays, and describe a sum rule estimate of the nonfactorizable 
_|- . contribution to the amplitude of B — * J/tpK. 
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1 Introduction 



During recent years a large amount of new experimental data on exclusive decays of charm 
and bottom hadrons has become available. In order to make optimal use of these data 
in determinations of fundamental parameters and tests of the standard model, one needs 
a quantitative understanding of the impact of strong interactions. This requires accurate 
calculations of hadronic matrix elements of weak operators in QCD beyond perturbation 
theory. Current approaches include lattice calculations, QCD sum rules, heavy quark ef- 
fective theory (HQET), chiral perturbation theory (CHPT), and phenomenological quark 
models. Each of these approaches has advantages and disadvantages. For example, quark 
models are easy to use and good for intuition. However, their relation to QCD is unclear. 
On the other hand, lattice calculations are rigorous from the point of view of QCD, but they 
suffer from lattice artifacts and uncertainties connected with the necessary extrapolations 
to the physical quark masses. Furthermore, effective theories are usually applicable only to 
a restricted class of problems, and sometimes require substantial corrections which cannot 
be calculated within the same framework. For example, HQET is very powerful in treating 
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b — > c transitions, but a priori less suitable for b — > u transitions, while CHPT is designed 
for processes involving soft pions and kaons. 

In this review, we focus on the method of QCD sum rules [|l], |2|, |3|. Proceeding from the 
firm basis of QCD perturbation theory this approach tries to incorporate nonperturbative 
elements of full QCD. Schematically, hadronic matrix elements of operators composed of 
quark and gluon fields are extracted from correlation functions of quark currents, rather than 
estimated directly in some models. One makes use of operator product expansion (OPE), the 
analyticity principle (dispersion relations), and S-matrix unitarity. In addition, one assumes 
the validity of quark-hadron duality in a rather strong form, sometimes called semi-local 
duality. The long-distance dynamics is parameterized in terms of vacuum condensates or 
light-cone wave functions. These new elements cannot yet be calculated rigorously in QCD, 
but are to be determined from experimental data in the one or other way. Nevertheless, 
because of the universal nature of the nonperturbative input, the sum rule approach retains 
its predictive power. Moreover, it is rather flexible and can therefore be applied to a large 
variety of problems in hadron physics. 

The application of sum rules to heavy quark physics has always been particularly success- 
ful. The more recent development of sum rule techniques for exclusive heavy meson decays 
again looks very promising. We shall substantiate this assessment for leptonic, semileptonic 
as well as nonleptonic decays of B and D mesons. Although the two flavours could be 
treated in parallel, we shall usually refer to B mesons, for definiteness. In most cases, it is 
obvious how to obtain the corresponding results for D mesons. 

The leptonic decay process B — > lui involves the simplest hadronic matrix element, that 
is the 5-meson decay constant We therefore have chosen the latter for our introductory 
study case. Furthermore, Jb is one of the important parameters of mixing and CP- violation 
in the B system, and a necessary input in order to eventually extract the CKM element V u b 
from future measurements of B — > fv T . Finally, fs also enters more complicated QCD sum 
rules for other hadronic properties of B mesons, e.g., form factors. 

The semileptonic decays B — > nlui and B — > plvi have already been observed experi- 
mentally. These exclusive modes are considered to provide interesting alternatives to the 
determination of V u b from inclusive b — > u transitions. However, truly competitive results 
can only be expected if the calculation of heavy-to-light form factors can be improved sig- 
nificantly. From the theoretical point of view, the B — > 7r and B —>■ p form factors are 
excellent examples for the application of light-cone sum rules. 

Exclusive nonleptonic decays are an even greater challenge to theory. Although over 
the years one has developed a qualitative understanding and even an amazingly consistent 
quantitative description of two-body decays 0, some features still lack a clear dynamical 
explanation. This concerns most of all the factorization of matrix elements of four-quark 
operators into products of matrix elements of quark currents. Since nonfactorizable am- 
plitudes are a priori expected to be channel- dependent, the apparent universality of the 
effective coefficients a% and ci2 constitutes a major puzzle. The decay mode B — > J/ipK 
containing two heavy quarks in the final state may play a special role in this respect and 
therefore shed some light on the factorization problem. We present this example also because 
it shows conceptual and technical limits of the sum rule approach. 

The content of this review is as follows. In section 2, we describe the derivation of 

2 A comprehensive study is presented by M. Neubert and B. Stech in this volume, ref. 0]. 
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the sum rule for fg from a two-point correlation function and summarize the analogous 
calculations for other B and D mesons. Also shown is a comparison with lattice results and 
experimental data. The straightforward generalization of the above method to three-point 
correlation functions is discussed in section 3 and applied to heavy-to-light form factors. In 
section 4, we explain an alternative approach based on light-cone expansion and employ it 
in section 5 to obtain sum rules for the B — > it form factors / + and f°. Recent results 
on the B — > p form factors are also mentioned here. Section 6 is devoted to sum rules for 
the strong B*Bn and D*Dir couplings. The estimates are used to normalize the single-pole 
approximation for / + and to calculate the width for D* — > Dn. In section 7, we present the 
sum rule predictions on the decay widths and distributions for B — ► nlu and B — > plv, and 
extract V u b from the comparison with data. Results are also shown on D — > irW. Section 
8 deals with the heavy mass expansion of the sum rules for heavy-to-light form factors 
and couplings, and with the asymptotic scaling behaviour. Finally, in section 9 we discuss 
the factorization hypothesis for B — > J / ipK and describe a first sum rule estimate of the 
nonfactorizable piece in the B — > J/ipK amplitude. 



2 Decay constants of B and D mesons 

Leptonic decays of B mesons are induced by the weak annihilation process bu — > lv[. The 
decay width 

T(B- - l-vfi = f^lKfclW^ f 1 -4) 2 fB (1) 
8tt V m B ) 

is suppressed by three small parameters: the CKM matrix element V u b, the lepton mass 
mi and the decay constant While the factor mf is enforced by helicity conservation, 
fs characterizes the size of the 5-meson wave function at the origin. Usually, the decay 
constant is defined by the matrix element of the relevant weak axial- vector current: 

(0 | q 1 , 1 ,b | B) = if B q, , (2) 

q carrying the flavour of the light quark constituent in the B, and q^ being the B four- 
momentum. Later we shall use the equivalent definition in terms of the matrix element of 
the corresponding pseudoscalar current: 

m b (0 | qi^ 5 b | B) = m%f B , (3) 

nib being the 6-quark mass. Because of the strong suppression of leptonic B decays only the 
mode B — ► fv T has a chance to be measured. However, this task will not be easy, even at 
future B factories. For this and other reasons which will become clear in the course of this 
review, accurate theoretical predictions of fs are indispensable. 

The QCD sum rule estimate of /# is based on an analysis of the two-point correlation 
function 1, 0|, § 



U(q 2 ) = % I d A xe iqx {Q | T{q(x)t-f 5 b(x),b(0)i-f 5 q(0)} | 0) . (4) 

By inserting a complete set of states with .B-meson quantum numbers between the currents 
in ([D one obtains the hadronic representation 

n( 2) = (Q\qil 5 b\B)(B\U l5 q\0) 
rn? B — q 2 
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„ < qt^ 5 b \h>< h \ h^ 5 q 0) 

+ 2^ 2 2 • ( 5 ) 

The decay constant fs appears in the first term contributed by the ground state S-meson. 
The second term takes into account the higher resonances and non-resonant states. The 
result ([5]) can be rewritten in the form of a dispersion relation: 



J irig s — q 



where the spectral density is given by 



, m %f 2 



p(s) = 8(s - m 2 B )^ + p h (s)&(s - 4) . (7) 



Obviously, the 5-function term on the r.h.s. of (|7j) represents the B meson, while p h (s) 
and Sq are the spectral density and threshold energy squared of the excited resonances and 
continuum states, respectively. In general, the dispersion integral ([]) is ultraviolet divergent. 
In order to make the integral finite, one subtracts a sufficient number of terms of the Taylor 
expansion of n(g 2 ) at q 2 = 0. In the case at hand, two subtractions are needed: 



iW) = rV) - n(o) - g 2 n'(o) = g 4 T . (8) 

Jm B S yS — q J 

Substitution of (|7|) in (||) yields 



n/(? 2 ) = ^ + I™ - n(o) - g 2 n'(o). (9) 



1 bIb 

ml(m 2 B — q 2 ) ' Js% s — q 2 
The subtraction terms can be removed by performing a Borel transformation with respect 



to q 2 



With 



and 



/ 2\ (n+1) / J \ n 

BM = Jz. ^^r- (s?) n(?2) s n(M2) ' (10) 



B M J^—\ — ^f'iV" \ en : 



( — \ P s/M 
^s-q 2 ) (k-l)\\M 2 J 



B AP (-q 2 ) k = (12) 

for k > one readily finds 

n(M 2 ) = n f (M 2 ) = H^ e -mi/M" + r p h< s ) e -s/M* ds (13) 

m£ Js% 

We see that Borel transformation removes arbitrary polynomials in q 2 and suppresses the 
contributions from excited and continuum states exponentially relative to the ground-state 
contribution. The second point is actually the main motivation for this transformation. 

At q 2 -C m 2 , below the poles and cuts associated with the resonances and continuum 
states, it is possible to calculate the correlation function IT(g 2 ) in terms of the quark and 
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gluon degrees of freedom in QCD. To this end, one expands the T-product of currents in 
(fD in a series of local operators Q d constructed from quark and gluon fields and normalized 
at the scale fi: 

% f d'xe^ T{q(x)ij 5 b(x), b(0)i l5 q(0)} = £ C d (q 2 , fi)Q d (fi) (14) 
J d 

with the sum running over the canonical dimensions d of the operators. The lowest- 
dimensional operators with d =0,3,4,5,6 are given by 

n d = 1, qq, G%G a ^, qaJ^GWq, {qT r q){qT s q) , (15) 

respectively. Here, A a are the usual Si7(3)-colour matrices, and is the gluon field 
strength tensor. T r denotes a certain combination of Lorentz and colour matrices. Note the 
absence of a colour-neutral, Lorentz-invariant operator with d = 2. 

The most important virtue of OPE is the possibility to separate long- and short- distance 
contributions to the correlation function 

Tl( q 2 ) = J2C d (q 2 ^)(0\nM |0> . (16) 

d 

While the strong interaction effects at momenta k 2 > fi 2 are included in the coefficients 
C d (q 2 ,fi), the effects at k 2 < fi 2 are absorbed into the matrix elements of the operators 
0<f(/jt). Thus, if \i 3> Aqcd the Wilson coefficients depend only on short- distance dynamics, 
and can therefore be calculated in perturbation theory, while the long-distance effects are 
taken into account by the vacuum averages (0 \ Q d \ 0) = (fl d ). These so-called condensates 
describe properties of the full QCD vacuum and are process- independent. At present, they 
can only be estimated in some crude approximations. For this reason, they are usually 
determined empirically by fitting sum rules to experimental data. It is essential for the 
whole approach that at q 2 -C mf the expansion ( |T6"D can be cut off after a few terms. The 
reason is that the higher the dimension of Q d , the more suppressed by inverse powers of 
m\ — q 2 is the corresponding Wilson coefficient C d . Most applications in the past include 
the set of operators given in ( p!5|) . 

The short-distance coefficients C d are calculated from the diagrams depicted in Fig. 1. 
Specifically, the coefficient Cq of the unit operator is obtained by contracting all quark fields 
in the correlation function (fD, and inserting the free quark propagators 

<0|TM*)g,(0)}|0> = iS^x) = 5 l3 J ^ c ^-*+^ , (17) 

where i,j are colour indices. For light quarks we put m q = 0. Diagrammatically, this 
approximation is represented by Fig. la. It constitutes the zeroth order approximation for 
the correlation function (^j) in QCD perturbation theory. The result can be written in the 
form of a dispersion relation: 

Tl Fig .aq 2 ) = C (q 2 ) = - r ds l ^$- (18) 

7T J ml S — q A 



4 




Figure 1: Diagrams determining the Wilson coefficients in the OPE of the two-point corre- 
lation function ($): C (a), C 3 (b), C 4 (c,d), C 5 (b,e), C 6 (b,f). Solid lines denote quarks, 
dashed lines gluons, wavy lines external currents. Crosses indicate vacuum fields. 



with 

ImC ( S ) = A (^11! . (19) 

87T S 



The subtraction procedure necessary to make the integral ( |T8| ) finite is postponed in antic- 
ipation of the Borel transformation fllPl). 

The coefficient C 3 of the qq operator in ( [Tfi| ) is obtained by treating the light-quark 
fields in ([|) as external vacuum fields and neglecting their momenta as compared to the 
momentum of the freely propagating off-shell b quark. The corresponding diagram is shown 
in Fig. lb. Substituting ( |17D in (f|) one gets 

IW fe (g 2 ) = / d'xe^iO | q(0)i«<l(Q)jP I OHTs^OrhsW , (20) 

where a,f3 are the spinor indices of the quark fields. Because of vanishing spin and colour 
of the vacuum, 

(0 | g(0) fa g(0) i/9 | 0) = ^^(qq). (21) 
Integration of ( pOl) over x then yields 

n Ftg . lb {<f)=Cz{q 2 )(qq) (22) 

with 

C 3 (q 2 )=m b /(q 2 -m 2 b ) . (23) 

The more complicated calculation of the coefficients C^^ from the diagrams of Fig. 1 is 
discussed, for example, in [HI. 




Figure 2: Feynman diagrams of the 0(a s ) correction to the Wilson coefficient Co(q 2 )- 



The accuracy of the OPE can be further improved by including perturbative QCD cor- 
rections to the Wilson coefficients corresponding to hard gluon exchanges in the diagrams 
of Fig. 1. Most important are the 0(a s ) effects on the coefficient Co 0, Q given by the 
two-loop diagrams of Fig. 2. Adding these corrections to the imaginary part of Co given in 
(|T§D, one has 

ImC.(.) = £^^(l + £/(.X)) , (24) 

where 



8ns \ 3ir 



f( S) ml) = - A +2U 2 (^)+\n—An 



4 s ml s — ml 

3, ml . s ml . s — ml ml . s 

+ - In + In ? + — b - In + In — = , (25) 

2 s — m£ s — m£ s mi s — mi m% 

\Ai{x) = — Jo — t), and m;, is the pole mass of the b quark. 

The complete result for the Borel-transformed correlation function (H) given in reads 

n«c D (M 2 ) = ^/ m , i^-j-Si- (l + ^/(«.™») «P ("^ 

+ ( - + i(^^G-) - ^ (l - ^) rnSfe) 

Wnosf m 2 b mj \ 2 \ ( m 2 h \ 

Here, use has been made of the conventional parametrization for the quark-gluon condensate 
density: 

(g V jG*" g ) = ^(w) , (27) 

and of vacuum saturation reducing four-quark condensates to squares of quark condensates 
with known coefficients [|TJ: 

(qT r qqT s q) = -^{(TrT r )(TrT s ) - Tr(T r T s )}(qq) 2 • (28) 

Now, by equating the hadronic representation of IT(M 2 ) given in ( |I~3D with the QCD 
result (|26|) one obtains an interesting sum rule which, however, does not yet allow to deter- 
mine Jb- This is due to the unknown spectral density p h (s) associated with the excited and 
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continuum states. In order to proceed, one assumes quark-hadron duality and substitutes, 
in (0), the perturbative spectral density: 

p h (s)e(s - 4) = -ImC {s)Q{s - s*) , (29) 

7T 

considering \[s^ as an effective threshold parameter of the order of the mass of the first 
excited B resonance. Subtraction of the integral in flnj|) from the corresponding term in 
Uqcd{M 2 ) then amounts to a simple change of the upper limit of integration, and yields 
the following sum rule for f B : 

3m? r s o , (s — ml) 2 /^ 4a, „, 9 A ( m% — s\ 



2 4 imt fa [s-mtY ( 4a, 2 \ 





M 2 



16n a s (qq) 2 ( m\ _ m\ \ \ ( m\-m\ \ 
'27 \ l 4M2 12M 4 / J y AP )' W 

The uncertainty due to the crude subtraction procedure should not be too harmful because 
of the suppression of excited and continuum states after Borel transformation. In principle, 
since the l.h.s. of ( p0[ ) is a measurable quantity, the scale dependence of the quark masses, 
the condensates and the running coupling on the r.h.s. must cancel. In practice, this can 
only be achieved approximately. 

The 6-quark mass appearing in the perturbative coefficient Co(g 2 , /i) is defined to be the 
pole mass, while the choice of rrib in the leading-order coefficients of the higher-dimensional 
terms in (|30|) is arbitrary. Following the usual procedure we take the pole mass 

m b = 4.7 ±0.1 GeV (31) 



everywhere in ([3y). The interval (pl|) covers the current estimates obtained from bottomo- 
nium sum rules ||. 

Turning to the condensates, it is important to note that the combinations rrib^qq) and 
(^ L G%G aflu ) are renormalization-group invariant. The numerical value of the quark con- 
densate density at /i = 0(1 GeV) is obtained from the PCAC relation [p], (|: 

mi GeV) = -^fe) =< -< 240 MeV » ' (32 > 



From that and the central value of ( |3~1~D converted to the running mass at 1 GeV, one finds 

m b (qq) ~ -0.084 GeV 4 . (33) 
The gluon condensate density is determined from charmonium sum rules [|TJ: 

<— G a uv G a H ~ 0.012 GeV 4 . (34) 
For the remaining parameters we take 

m 2 (l GeV) ~ 0.8 GeV 2 , (35) 



7 



as extracted from sum rules for light baryons 
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a s (qq) 2 = 8- 10" 5 GeV 6 , (36) 

as given in 0. The scale dependence of m^qq) and a s {qq) 2 is negligible. The numerical 
uncertainties on the condensates vary from 10 % to about 50 % or more. However, they 
influence the final result on rather little. 

This leaves us with the choice of scale in the running coupling a s . As the average 
virtuality of the quarks in the correlator (|j) is characterized by the Borel parameter M 2 , it 
is reasonable to take 

= 0(M 2 ) . (37) 

Finally, we have to determine the values of M 2 for which the sum rule ( P0| ) can be trusted. 
On the one hand, M 2 has to be small enough such that the contributions from excited and 
continuum states are exponentially damped and the approximation ( p9|) is satisfactory. On 
the other hand, the scale M 2 must be large enough such that higher- dimensional operators 
are suppressed and the OPE converges sufficiently fast. Consequently, there is at best a 
finite range of M 2 in which a given sum rule is valid. Moreover, in this range the numerical 
result should be stable under variations of M 2 . Whether or not these requirements can be 
met has to be investigated for each case separately. For the sum rule ([30]) we have checked 
that with 

3 GeV 2 < M 2 < 5.5 GeV 2 (38) 

the necessary conditions are fulfilled: the quark-gluon condensate contributes less than 
15 %, the contributions from the gluon condensate and from the four-quark operators are 
negligibly small, and the excited states and continuum contribute less than 30 %. Moreover, 
the change of Jb due to variation of M 2 in the range fl38|) is indeed very small. 

The value of Jb derived from the sum rule (|30|) is given in Table 1. The value changes 
from 210 MeV to 150 MeV when is varied in the range fl3T|) from 4.6 to 4.8 GeV, and 
at the same time from 37 to 33 GeV 2 . The correlation between and improves 
the stability of the sum rule under variation of M 2 . As compared to the above uncertainty, 
the uncertainties from the condensates are negligible. Note that the effect of the 0(a s ) 
correction (^4j) to Co is sizeable. Without this correction the value for fs is much lower: 

J B = f B (a s = 0) = 140 ± 30 MeV . (39) 

Also shown in Table 1 is the corresponding prediction on fo- It follows from (|30|) after 
replacing by the pole mass m c = 1.3 ± 0.1 GeV of the c quark, and s B by the effective 
threshold = 6 1 GeV 2 in the D-meson channel. The allowed interval in the Borel mass 
is found to be 1 GeV 2 < M 2 < 2 GeV 2 . For later use, we again quote the result without 
the 0(a s ) correction: 

J D ee f D (a s = 0) = 170 ± 20 MeV. (40) 

In addition, Table 1 allows comparing QCD sum rule estimates with lattice results and 
with the available experimental data. The mutual agreement, within the uncertainties of 
the two theoretical methods, is satisfactory. On the other hand, the data are just beginning 
to challenge theory. 



Table 1: Decay constants of B and D mesons in MeV. 



Method 


Ref. 


Ib 


/-B 


Id 


Id 


QCD sum rules 


re 
] 


th 
;vi 

1 


is 

ew 

a) 


180 ± 30 
175 


210 


190 ± 20 
180 ± 10 


220 ± 10 


Lattice 




11. 

E 
i 


i 

5] 


175± 25 
180± 32 


200± 25 
205± 35 
196li 


205± 15 
221± 17 
19113 


235 ± 15 
237± 16 
2061^ 


Experiment 


[] 


ii 

L7 


3] 








< 310 
241 ± 21± 30 



a ) update of results of [g, 0, |TJ taking m b = 4.67 GeV, s% = 35 GeV 2 , 
and m c = 1.3 GeV, = 5.5 GeV 2 . 

b ) world average 



It should be mentioned that the S-meson decay constant has been derived also from a 
two-point sum rule in the heavy quark limit |18| using HQET and including 1/rriQ corrections 



[|T^, PC|| . m this approach, the scale in a s is lower than (|3?D, of the order of the reduced 
mass in the 5-meson bound state. This leads to a somewhat larger value for close to 
the upper end of the range obtained from the sum rule Q3T]|). For more details, one may 
consult the original papers quoted above. 

Finally, for the calculation of the B*Bir and D*Dir couplings in section 6, we will need 
the B* and D* decay constants. Since they can be estimated from similar two-point sum 
rules as Jb and fc, we briefly summarize the result here. The decay constant of the B* is 
defined by the matrix element 

(0 | qinb | B*) = vriB*fB*t l i , (41) 

e M being the B* polarization vector. From the correlation function of the vector currents 
(Pipb and b-f^q similar to (f|), and with the same approximations as in (|30"D, one derives the 
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following sum rule: 



A-4. - i £f ( 2 + * ) (l + £ r ( , exp ft-) 

+ { - (l - ^) + C 4 (M 2 )(^G^> 

+ C 6 (M 2 K(gg> 2 } exp (^^) • (42) 

The explicit expressions for the coefficient /* of the 0(a s ) correction, as well as the coeffi- 
cients C4 and Cq of the numerically insignificant contributions from the gluon and four-quark 
condensates are given in ||. The corresponding sum rule for f&* can be directly inferred from 
( |42|) . For the parameters one may take the same values as in the case of the pseudoscalars. 
For later use we give below the leading-order results in which the 0(a s ) loop-corrections to 
the coefficient of the unity operator are not included: 

f B » = f B , (a s = 0) = 160 ± 30 MeV , (43) 
Jd* = fo* (as = 0) = 240 ± 20 MeV . (44) 



3 Three-point sum rules for heavy meson form factors 

We next turn to the more complicated matrix elements of hadronic transitions induced by 
weak currents. An important example is the B — ► 7r transition. The relevant matrix element 
is parametrized by two independent form factors: 

(n(q)\u^b\B{p + q)) = 2f+(p 2 )q, + (/ V) + T (/)) P, (45) 

with p + q, q and p being the B and ir momenta, and the momentum transfer, respectively. 

Generalizing the method employed in section 2 one starts from the following vacuum 
correlation function of three currents: 

T a ,(p,q) = -J rf 4 xAe^ + ^(0|T{%) 7a 75^(y)^(x)7^(^)J(0)^(0)}|0) , (46) 

where the B and 7r meson states of ( [l"5|) are replaced by the corresponding generating 
currents. Insertion of the complete sets of hadronic states carrying B and it quantum 
numbers, respectively, yields the double dispersion relation 

Ta (p q ) = m lfBf*Q« (V + (P 2 )% + (f + (p 2 ) + f-(p 2 ))p») 



m b (m B - (p + q) 2 )(ml - q 2 ) 

Pa»( S U S 2,P 2 ) 



+ 



/ / ds 1 ds 2 -. — 

J JSia [Si 



(p + q) 2 )(s 2 -q 2 ) 

+Pitf,P>) I d Sl ph ^ Sl fl + P 2 ((p + q) 2 ,p 2 ) I ds / 2 ^ S2 { ] . (47) 
JSi si — [p + qy Jt.2 s 2 - q z 
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In the above, the matrix elements (Q), fl45|) , and 

(0 | dlal^u | 7r(g)) = if v q a (48) 

have been used. While the first term on the r.h.s. of (|47|) is the the B and 7r ground state 
contribution, the integral over the double spectral density p^ takes into account the con- 
tributions of higher resonances and continuum states in the B and 7r channels. S 12 denotes 
the integration region in the (si,S2) plane. B-functions defining the actual thresholds are 
implicitly contained in p\ . The additional single dispersion integrals multiplied by polyno- 
mials P\ t 2 arise from subtractions, similar to the polynomial in (^). They vanish after Borel 
transformation. 




Figure 3: Diagrams determining the Wilson coefficients of the OPE of the three-point cor- 
relation function The symbols are as in Fig. 1. 

At (p + q) 2 <C ml, g 2 < 0, and p 2 < m\ — ZmbXi X being a m&-independent scale of 
order Aqcd, the correlation function T afJi (p,q) can be approximated by the first few terms 
in the OPE of the T-product of currents in (fJ6f), in analogy to ([16]): 

T^(p,q) = J2Cr(p,q,^M). (49) 
d 

The restriction on p 2 is necessary in order to stay sufficiently far away from the physical 
states in the bu channel, most notably the B* . For d < 6 the short- distance coefficients 
can be calculated from the diagrams shown in Fig. 3. The calculation follows the 
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procedure outlined in the previous section. However, the kinematics is more complicated 
because of the presence of two independent external four-momenta. The explicit expressions 
can be found in pl| |23|. Decomposing as well as the spectral densities p\ in the 
independent tensor structures, 



CT (p, q) = C d ((p + q)\ q 2 ,p 2 )q a (2q + pf + - 

p^(si,s 2 ,p 2 ) = p h (si } s 2} p 2 )q a (2q + p)^ + ... . 
and equating the corresponding invariant functions in (|47D and (fl9|), 

™ 2 B fBfJ + (p 2 ) 



(50) 
(51) 



one obtains 



+ 



m b (m 2 B - (p + q) 2 ) (m 2 - q 2 ) 'I Is 12 dsids2 ( Sl - (p + q) 2 ) (s 2 - g 2% 

= £c d ((p + g)^v^)<^(p)> 



p h {s ll s 2l p 2 ) 



+ ... 



(52) 



where the ellipses denote the subtraction terms. Since light quarks are taken massless, we 
put m n = for consistency. 

Analogously to fl29"l) , the double spectral density p h is approximated by 



P h (si, s 2 ,p 2 ) = — Im sl Im S2 C (si, s 2 ,p 2 )Q(s 1 - Sq)Q(s 2 - sj) , 



7T 



(53) 



where Co is the perturbative coefficient calculated from Fig. 3a, and Sq is the effective 
threshold parameter in the 7r-meson channel. After Borel transformation (|I]) with respect to 
the variables (p+q) 2 and q 2 the subtraction terms disappear from (|52|) and the contributions 
from excited and continuum states become exponentially suppressed. The final sum rule for 
/ + takes the following schematical form [pi , p^ , [23| , 24]: 



f + iP 2 



m 



B 



m 2 B f B U eXP \m} + M 2 



7T- 



Sl 



s 2 



Sl2(-ff,»J) 



dsids 2 Im Sl Im S2 Co(si,S2,P ,p)exp -770 _ 772 



M 2 



+ ^C d (M 1 2 ,M 2 2 ,p 2 ,p)(^(p)) 

d=3 



(54) 



with M 2 and M| being the Borel parameters in the B and tt channels, respectively. The 
appearance of the threshold parameters s$ and Sq reflects the continuum subtraction. Sim- 
ilarly as for the decay constants, there may be important perturbative QCD corrections 
to the short-distance coefficients, in particular, to Co. To our knowledge, the two-loop 
diagrams corresponding to hard-gluon exchanges in Fig. 3a have not yet been calculated. 
The numerical estimate for f + {p 2 ) obtained from Q53D is plotted in Fig. 4 in com- 



parison with other predictions. The values for m^, Sq, and the vacuum condensates are 
already stated in the previous section; for fs the leading-order estimate (B^) is used for 
consistency. Furthermore, the pion decay constant is f n = 132 MeV, while the threshold 
parameter Sq = 0.75 GeV 2 is inferred from the two-point sum rule for the correlation 
function (0 | T{d , y a 'y 5 u(x), u^^d^O)} | 0). 
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Similar three-point sum rules have been applied to a variety of hadronic problems in 

I 



eluding the pion form factor |27 
decays |3C 



2j§, and rare B 



radiative charmonium transitions 
32] . However, there are theoretical difficulties which put some doubts on the 
Most disturbing in the case of heavy-to-light form factors is the 

More specifically, the coefficients of the 



reliability of this method. 

breakdown of the OPE in the heavy mass limit, 
subleading quark and quark-gluon condensate terms grow faster with mj than the coefficient 
of the leading unity operator, that is the perturbative contribution. In contrast, the heavy 
mass limit of the two-point sum rules for decay constants discussed in the previous section is 
completely well-behaved. A more detailed discussion on this issue can be found in |33"| , |34| . 

Although three-point sum rules still remain a useful calculational tool for selected prob- 
lems, the above difficulties call for more consistent and reliable sum rule methods. One 
important development is explained in the next section. 



4 Light-cone expansion 

For processes where a light meson such as a tt, K, or p is involved, there is an interest- 
ing alternative to the short-distance OPE of vacuum-vacuum correlation functions in terms 
of condensates, namely the expansion of vacuum-meson correlators near the light-cone in 



terms of meson wave functions [36|, [37]. The latter are functions of the light-cone momen- 
tum fractions carried by the constituents of a given meson. Similarly as the deep-inelastic 
structure functions, the wave functions can be classified by the twist of the corresponding 
composite operators. While the light mesons are taken on mass shell, the heavy meson 
channels are treated in the usual way: choice of a generating current, contraction of the 
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heavy quark fields, dispersion relation, Borel transformation and continuum subtraction. 



The light-cone variant of QCD sum rules suggested in |3S, |39|, |4(J allows to incorporate 
additional information about the Euclidean asymptotics of correlation functions in QCD for 
arbitrary external momenta. Moreover, it avoids the problems of the three-point sum rules 
mentioned at the end of the previous section. 

Here, we explain the general idea using the B — > n transition element (45) as an example. 
The detailed derivation of the light-cone sum rules for / + and /~ [41, 42, H3| is postponed 



to the next section. The starting point is the vacuum-pion correlation function 
F,(p,q) = % f d 4 xe ipx (n(q) | T{u(x) lti b(x), b(0)i l5 d(0)} | 0) 



(55) 



F(p 2 , (p + qf)q^ + F{p\ (p + qf)p^ . 



Since the pion is on-shell, q 2 = m 2 vanishes in the chiral limit adopted throughout this 
discussion. For the momenta in the bd and ub channels we respectively require (p+q) 2 m 2 
and p 2 < m 2 — 2mf } x, as before. Contracting the 6-quark fields in (|55|) and inserting the free 
6-quark propagator (|17|) , one gets 



d 4 x d 4 k 



(2vr) 4 (m2 - k 2 ) 
+ r(7T(g)|n(x) 7M 7,75rf(0)|0)) 



i( *- k > (m b (n(q)\u(x)^d(0)\0) 



(56) 



This contribution is diagrammatically depicted in Fig. 5a. 




Figure 5: Diagrammatic representation of the correlation function (\53j): terms involving 
two-particle (a) and three-particle (b) wave functions. Solid lines represent quarks, dashed 
lines gluons, wavy lines external currents, and ovals light-cone wave functions of the pion. 

Let us first focus on the matrix element 

(7r(g)|n(x)7 M 7 5 rf(0)|0), (57) 
and expand the bilocal quark-antiquark operator around x = 0: 

u(xh^ 5 d(0) = ]T ^w(0)(B -xY^diO) . (58) 
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The matrix elements of the local operators can be written in the form 



(n(q)\u D ai D a2 ... D ar 7^75^)0) = (i) r q ti q ai q a2 --.q ar M r + ... , (59) 

where the ellipses stand for additional terms containing various combinations of the metric 
tensor g ai0lk . Note that the lowest twist in (|59|) is equal to 2. Substituting (|58|) and (|59|) 



in the first term of ([x]), integrating over x and k, and comparing the result with (|55j), one 
obtains 

F(p 2 , (p + qf) = i^~ 2 £ f M r (60) 

with 

2(p-g) (p + g) 2 -p 2 , , 

£= — ~ = 2 2 — ■ ^ 61 ) 

mi — p z mi — p z 

Now, one immediately encounters the following problem. If the ratio £ is finite one must keep 
an infinite series of matrix elements in (^0]) . All of them give contributions of the same order 
in the heavy quark propagator l/{m 2 —p 2 ), differing only by powers of the dimensionless 
parameter £. In other words, the expansion ([58]) is useful only in the case £ — > 0, i.e., when 
p 2 ~ (p + q) 2 , or equivalently when the momentum of the light meson vanishes. Under this 
condition, the series in (|5Up can be truncated after a few terms involving only a manageable 
number of unknown matrix elements M r . However, generally, when p 2 ^ (p + q) 2 one has 
to sum up the infinite series of matrix elements of local operators in some way. 

One possible solution is provided by expanding the bilocal operator (|58D near the light- 
cone. For the matrix element (Elt) the first term of this expansion is given by 



(7r(q)\u(x)^j 5 d(0)\0) x2=0 = -iq^ J due luqx <p w (u) . (62) 

The function ^(w) is known as the twist 2 light-cone wave function of the pion. It is 
normalized to unity, as can be seen by putting x = in (|62|). As already mentioned, (p^ 
represents the distribution in the fraction of the light-cone momentum go + 93 of the pion 
carried by a constituent quark. The path-ordered gluon operator 

Pexp ^ig s J da x IM A fi (ax)^ (63) 

necessary for gauge invariance of the matrix element (0), is unity in the light-cone gauge, 
x^A^ = 0, assumed here. Substitution of (^) in fl5"6|) and integration over x and k yields 

iV, (p + q) 2 ) = m b U I' 2 dU ^ ] . (64) 

Jo ml — [p + uqy 

We see that the infinite series of matrix elements of local operators encountered before in 
(|60"D is effectively replaced by a wave function. If one expands (B4|) in q, 



F(p 2 , (p + q) 2 ) = m b U £ , ( ?' q L +1 I" du u r ^(u) , (65) 



3 Twist is defined as the difference between the canonicai dimension and the spin of a traceless and 
totally symmetric local operator. 
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and compares the above expression with (|60"D , one can directly read off the relation between 
the matrix elements M r defined in (|59"D and the moments of (p n (u): 



M r 



i/ w / du u r (p^(u). 



(66) 



It should be noted that the pion wave function ip n (u) is a universal function. It is the 
same quantity which enters, for example, the 7r°7*7* form factor ||35|| . This universality is 



essential for the whole approach, similarly as the universality of the vacuum condensates 
for the short- distance sum rules. 

Including the next-to- leading terms in x 2 , the light-cone expansion of the matrix element 
reads 



(7r(q)\u(x)^ 5 d(0)\0) = -iq^ / due** UpM + x 2 9l (u] 

J 

+fJx fl -^) ^ due^ x g 2 (u), 
\ qx I Jo 



(67) 



where g\ and g 2 are twist 4 wave functions. Proceeding to the second term in ( p6|) and using 
the relation 



one encounters two further matrix elements: 



(6* 



and 



(n(q) | u(x)i-f 5 d(0) | 0) = / du e iuqx <p p (u) 

Jo 



(vr(g) | u(x)a^ 5 d(0) | 0) = i(q^x u - q^x^)^ 1 - [ du e Mqx (p a (u) 



6 Jo 



(69) 



(70) 



with Utt = m 2 /(m u + rrid). Only the leading terms in the expansion are considered here. 
They have twist 3 and are parameterized by the wave functions (p p and (p a . 



Beyond twist 2, one should also take into account the higher-order terms j44| resulting 
from the contraction of the 6-quark fields in the correlator (|55): 



(0|T{6(x)6(0)}|0> 



iS b (x) - ig s 



d 4 k r 1 

' 1 dv 



(2tt 



1 ft + m b 



+ 



1 



mr — k'- 



-.vXfjG^^x)^ 



(71) 



where S is the free propagator given in (|ITD , and G pu = G c pv — with tr(A a A ) = 25 . Inser- 
tion of the gluonic term in the correlation function ([55]) yields the contribution represented 
by the diagram Fig. 5b: 



FP(p,q) = ig s 



d k d x dv 
(2vr) 4 (m2 - k 2 ) 



^ p ~ k)x (7i\u(x)jJvx p G p \vx)-f X 



+ » 



1 ft + m b ri>x 

2 ml — k 2 



G px (vx)a p x )t5^(0)|0) 



(72) 
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Use of the identities 



(73) 



and 



lulvOpX = {op\9v P - (ypp9u\) +i(gp\9up- QppQux) 



(74) 



leads to the following matrix elements of quark-antiquark-gluon operators which bring three- 
particle wave functions into the game: 

(ii\u(x)a flu j 5 g s G a p(vx)d(0)\0) = 
= ifa7r[(.QaQii9pv - Wp.9olv) ~ (q a qugpfi - W^ap)] J Va t (f 3n (ai)e iqx( - ai+vas) , (75) 



(ir\u(x)j IJ ,j 5 g 8 G a p(vx)d(0)\0) 

■EaQp \ 



qp y9ap 



qx 



qx / 



+ fit — (<laXf3 ~ qpXa) \ T>Q>i ip\\(oii)e 



qx 



Jqx(a 1 +va 3 ) 



(76) 



(w\u(x)"/ li g s G a p(vx)d(0)\0) 



q^ 



qx 



q a g?p 



qx / 



Von ip±{ai)e 



iqx(ot\+va3) 



+ if^{q a xp - qpx a ) \ Va % ^(a 4 )e^ ai+ ^ 
qx 



(77) 



with G a/ 3 = \e a p aT G (TT , and Vai = da\da2da 3 5{l — a.\ — ai — a 3 ). The wave function 
y37r(tti) = <P3n(ai, «2, as) has twist 3, while ip±, <p\\, (p±_ and (p\\ are all of twist 4. Gluons 
emitted from the light quark lines in Fig. 5a are effectively taken care of by the twist 3 and 



4 wave functions as was shown in |39], [45], [±6|. Components of the pion wave function with 
two extra gluons, or with an additional qq pair are neglected. 

Finally, a comment is in order concerning the perturbative corrections to the correlation 
function (p5]) from hard gluon exchanges. The first-order diagrams are shown in Fig. 6. 



Returning to the twist 2 approximation ( 64]) of the invariant function F, one can more 
generally write it as a convolution of a hard scattering amplitude T with the twist 2 wave 
function: 



F(p ,{p + q) ) = U / dwp w (u, n)T(p ,{p + q) ,u, fi) 
Jo 



(7f 



where T is given by 



T(p 2 , (p + q)\ u, ft) = T (p 2 , (p + q) 2 , u) + ^Up 2 , (p + q) 2 , u, ft) + 0(a 2 s ) . (79) 
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Figure 6: Diagrams of the 0(a s ) corrections to the correlation function (j5. 



As well known, the renormalization procedure and the factorization of the collinear loga- 
rithms generated by gluon radiation induce scale dependences. In ( |T8| ) we have chosen a 
common scale /x, for simplicity. 

In leading-order approximation (LO), the lowest-order amplitude 



T (p ,{p + q) ,u) 



m l — P 2 (^ — u) — (p + o) 2 u 



0) 



is convoluted with the wave function obeying the Brodsky-Lepage evolution equation |35 

i 

d(p w (u, /i) / din n = J dwV(u,w)(f 7T (w, fx) (81 
o 

with 

a s (/i)C F 



V(u, w) 
and 



7T 



i + e(u - w) + - i + 

1 — w \ u — wj w V w — u 



Q(w — u) 



R(u, w) 



R(u, w) — 8{u — w) J R(v, w)dv 



82 



(83) 
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This equation effectively sums up the leading logarithms to all orders. It can be solved by 
expanding (p n (u, p) in terms of Gegenbauer polynomials in which case the coefficients are 
multiplicatively renormalizable [|35|j . 

The next-to-leading order approximation (NLO) is given by the hard amplitude (|7"9| ) 
including the first-order term T\ f|7|, 48|, together with the wave function <p n (u, p) solving 



the evolution equation in next-to- leading order p9j . To date, the NLO program is completed 
only for the twist 2 contribution to the invariant amplitude F. The relevant wave function 
parameters and correction terms are given in Appendices 1 and 2. 



5 Light-cone sum rules for heavy- to- light form factors 

In order to determine the B — > n form factors f + {p 2 ) and f~(p 2 ) from the correlation 
function (|55|) , we employ a QCD sum rule with respect to the 5-meson channel following 
essentially the same steps as in the derivation of the sum rule for f B in section 2. The 
hadronic representation of (|5f|) is obtained by inserting a complete set of intermediate states 
with 5-meson quantum numbers: 

(n\u %l b\B)(B\U l5 d\0) 
m%-(p + q) 2 

y- (?r I u-f^b \h)(h\ bi~f 5 d | 0) 
Using the matrix elements (pi) and (|45|), and representing the sum over excited and contin- 



uum states by a dispersion integral with the spectral density p 1 ^ = p h q fl + p h p^, one obtains 
the following relations for the invariant amplitudes F and F: 

2mlf B f+(p 2 ) | J p h (p 2 ,s)ds 
m b (m 2 B - (p + q) 2 ) J s - (p + q) 2 



m {P + ,)•) = nteWhim + J^^L . ( sa) 

m b {m B - (p + qY) J s - (p + q) z 



Similarly as in (p5|) , the integrals over p h and p h are again approximated by integrals over 
the corresponding spectral densities calculated from the light-cone expansion of (55) using 

p h (p 2 , s)Q(s - s h ) = -lmF QCD (p 2 , s)Q(s - s s ) , (87) 

7T 

and the analogous relation for ph and IirFqcd- The calculation of Fqcd and Fqcd has 
already been outlined in the preceding section. With (87) it is now straightforward to sub- 
tract the contribution of the excited and continuum states from the corresponding integrals 
on the l.h.s. of (|85f) and (|86|). After performing the Borel transformation in (p + q) 2 , one 
finally arrives at the sum rules 



s 
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fB(f + (p 2 )+f-( P 2 )) 



m b 



nm 



B 



ImF Q cD(p 2 ,s)exp 



m B — s 
M 2 



ds . 



39) 



From the above, the form factor / follows directly by subtraction, while the scalar form 
factor f° is given by 



/V) = f + ( P 2 ) + 



p 



m 2 B — m 2 



f-(p 2 



(90) 



The remaining task is to complete the calculation of Fqcd and Fqcd, and to determine 
the imaginary parts. We proceed along the lines explained in section 4. The contribution 
of two-particle wave functions is derived from the diagram Fig. 5a inserting the matrix 
elements (p7D, (1691), and (fTOh in (pra). This results in 



f$ d (p 2 ,(p + q) 2 ) = U 



du 



Vp( ) + g ^ + m 2 _ ^ _|_ 



m b ~~ (p + u q¥ 

p 2 + ml 



+m b 



2ug 2 (u) 



m b ~ (P + u l) 2 ( m b ~ ip + uq) 2 )' 



dvg 2 {v) 



(91) 



^(P 2 ,(p + g) 2 ) = /ir 



mf — (p + uq) 2 



6u 



mt 



p 



m b ~ i.P + M '?) 2 



2m b g 2 {u) 



mt 



(p + uq) 2 



Furthermore, the three-particle contribution is obtained from Fig. 5b using 
matrix elements ( [75]) to (f77|): 



^(P 2 ,(P + ?) 2 ) 



(if / Pa 



4/3 7 r^37r(ai) V (P'?) 



,(fc) 



k [to 2 - (p + («i + wa 3 )g) 2 ] 2 

2^±(oii) - f^joj) + 2<p ± (ai) - (f>\\{OLi 
[ml - (p + (ax + v a 3 )q) 2 ] 2 

. 



(92) 
and the 



(93) 

Q CD (P 2 ,(P + <1) 2 ) = 0. (94) 

It is interesting to note that there are no contributions from twist 2 and three-particle wave 
functions to Fq C d, and hence to f + + f~- Including all operators up to twist 4 one has in 
total 

F Q cd(p 2 , (p + qf) = F$ D (p 2 , (p + qf) + F^ CD (p\ (p + q) 2 ) , (95) 



F^ n {p\{p + qf). 



Fqcd( P 2 , (P + q) 2 ) = F$UP\ (P + lY)- (96) 

We see that every power of x 2 in the light-cone expansion of the integrand in (]55j) leads to 
an additional power of the denominator m 2 — (p + uq) 2 in Fqcd and Fqcd- This justifies 
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the neglect of higher-twist operators, provided (p + q) 2 <C m 2 and p 2 < ml by a few GeV 2 . 
More definitely, numerical studies indicate that contributions beyond twist 4 can safely be 
neglected. 

The imaginary parts of the two-particle functions Fqq D and Fqcd are relatively easy 
to find. Terms proportional to l/[ml — (p + uq) 2 ] can be directly converted into disperson 
integrals with respect to (p + q) 2 by changing the variable u into s = (ml — p 2 )/u + p 2 . 



Subtraction of the continuum simply shifts the lower limit of intergration in (pi]) and in (|92|) 
from to A = {ml—p 2 )/(sQ—p 2 ). Terms proportional to higher powers of l/[ml — (p+uq) 2 } 
need to be partially integrated. After continuum subtraction, this leads to additional surface 
terms which are not completely negligible. A similar procedure is possible in the case of 
the three-particle function Fq^ d . However, since this contribution is anyway only a small 
correction, we neglect the corresponding surface term, for simplicity. 

In summary, fl5%|), (|9"T|), (|9"3"|), and (|9l)| ) yield the following sum rule for / + , or more 
precisely, for the product fsf~ 



41, 42 



fsf + (p 2 ) 



firml 
2m| 



exp 




ml — p 2 {l — u) 
^M 2 



x (p n (u) H 

V m b 



u(p p (u) + 



Am 2 bgi (u) 2 



u 2 M 4 



uM 2 



<fa(u) 



g 2 {v)dv 1 + 



1 + 



9 9 s 

m^ + p 



m 2 + p 2 



./;;(/r..\/ 2 ) • 




2uM 2 j 
+ t+(s*,p 2 ,M 2 ) 



(97) 



Here, t + denotes the surface term just mentioned, Jq the contribution from the three- 
particle wave functions, and 6 + the 0(a s ) correction to the leading twist 2 term [[47], |4~8f 
Explicit expressions can be found in Appendix 2. Similarly, the sum rule for / + + /" 
is obtained from (^) and d96|): 



/ B (/ + (/) + r(p 2 )) 



ni 



exp 




exp 



mt 



p (1 - -u) 



wM 2 



1 - 



+ 



2m b g 2 {u) 



M 



(9^ 



6m V" mM 2 / /i^M 2 

with the surface term t 1 * 1 being given in Appendix 2. Perturbative QCD corrections to this 
sum rule are still missing. 

In the following, we shall illustrate the above results numerically. For this purpose, 
it suffices to work in LO. Strictly speaking, in the absence of a complete NLO analysis 
this is also more consistent. Furthermore, for nib and J'b we choose the values given in 
( |3T|) and (p9|), respectively, for the threshold parameter we take Sq = 35 ± 2 GeV 2 . The 
pion wave functions are listed in Appendix 1. With this input we have checked that for 
M 2 = 10 ±2 GeV 2 the twist 4 corrections are very small, and the contributions from excited 
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4 8 12 16 

p 2 [GeV 2 } 

Figure 7: B — > 7r /orm factors obtained from light-cone sum rules in LO. 

and continuum states do not exceed 30%. Up to p 2 ~ 18 GeV 2 , the sum rules are also quite 
stable with respect to a variation of M 2 . However, at larger momentum transfer the sum 
rules become unstable, and the twist 4 contribution grows rapidly. This is not surprising 
because the light-cone expansion and the sum rule method are expected to break down as 
p 2 approaches m 2 . If not stated otherwise, the numerical results given below are obtained 
from the central values of the input parameters, which we refer to as our nominal choice. 

The momentum dependence of the form factors / + and / + + /~ can be seen in Fig. 7. 
Specifically, at zero momentum transfer we predict 

/ + (0) = 0.30, 

/+(0) + /-(0) = 0.06. (99) 

If the known 0(a s ) corrections are included in the sum rules for /s/ + (p 2 ) and fs, one gets 
instead 

/ + (0) = 0.27. (100) 

Of course, a most important question concerns the reliability of these predictions. Below 
we comment on the main sources of uncertainties. Allowing the parameters to vary within 
the ranges given above and in Appendix 1 we use the deviations of / + from the value 
obtained with the nominal choice of parameters as an uncertainty estimate. 

(a) Borel mass parameter 

The dependence of / + on M 2 is illustrated in Fig. 8a. In the allowed interval of M 2 , f + 
varies by only ±(3 to 5) %, depending on p 2 . 

(b) 6-quark mass and continuum threshold 

Fig. 8b shows the variation of / + with rrif, keeping all other parameters except fixed. The 
value of /b is taken from the sum rule (|30"D dropping the 0(a s ) corrections. The analogous 
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Figure 8: Variation of the sum rule prediction on f + (p 2 ) with the Borel parameter (a), the 
b-quark mass (b), and the threshold parameter Sq (c). Considered are three typical values 
of momentum transfer: p 2 = (solid), p 2 = 8 GeV 2 (long -dashed), and p 2 = 16 GeV 2 
(short-dashed). 

test for Sq is performed in Fig. 8c. If Sq and m b are varied simultaneously such that one 
achieves maximum stability of the sum rule for f B , the change in / + is negligible at small 
p 2 rising to about ± 3 % at large p 2 . 

(c) higher-twist contributions 

No reliable estimates exist for wave functions beyond twist 4. Therefore, we use the magni- 
tude of the twist 4 contribution to / + as an indicator for the uncertainty due to the neglect 
higher-twist terms. From Fig. 9a we see that the impact of the twist 4 components is 
comfortably small, less than 2 % at low p 2 and about 5 % at large p 2 . This suggests a 
conservative estimate of ± 5 % due to unknown higher- twist. 

(d) light-cone wave functions 

The asymptotic wave functions and the scale dependence of the nonasymptotic coefficients 
are given in perturbative QCD. However, the values of these coefficients at a certain input 
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Figure 9: B — > 7r /orm factor f + : (a) individual contributions of twist 2 (dashed), 3 (dash- 
dotted), and 4 (dotted), and the total sum (solid); (b) wave functions with (solid) and without 
(dashed) nonasymptotic corrections. 
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scale no have to be determined empirically. They are presently only known with considerable 
uncertainties. In order to clarify the sensitivity of / + to nonasymptotic effects, we put the 
latter to zero and compare the result in Fig. 9b with the nominal prediction. The change 
amounts to about -10 % at small p 2 and +10 % at large p 2 , while the intermediate region 
around p 2 ~ 10 GeV 2 is very little affected. Since this exercise is rather extreme, a more 
careful estimate should give a smaller number. Note that the sum rule being dominated by 
convolutions of smooth amplitudes with normalized wave functions, is actually expected to 
be relatively insensitive to the precise shape of the latter. 
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Figure 10: B — > 7r form factor f + (p 2 ) in leading twist 2 approximation: LO (dashed) in 
comparison to NLO (solid). 



(e) perturbative corrections 
The 0(a s ) corrections to the leading twist 2 term in the sum rule (p7[) for fsf + and to 
the sum rule ( |30"|) for both being about 30 %, cancel in the ratio. The net effect is the 
unusually small correction to / + shown in Fig. 10. This result H7L HS1 eliminates one main 



uncertainty. Unfortunately, perturbative corrections to the higher-twist terms have not yet 
been calculated. Considering that the twist 3 terms contribute about 50 % to the sum rule 
for and assuming again radiative corrections of about 30 % but no cancellation, one 

has to face a remaining uncertainty of about 15 %. This, however, may be an overestimate. 
More importantly, with some major effort this deficiency can be cured. 

In summary, the present total uncertainty in f + (p 2 ) is estimated to be about 25 %, if the 
uncertainties (a) to (e) with the exception of (d) are added up linearly, and about 17 % if 
they are added in quadrature as is often done in the literature. Once the perturbative QCD 
correction to the twist 3 term is calculated, this uncertainty, which mainly concerns the 
normalization, can be reduced to 10 %. In addition, there is a shape-dependent uncertainty 
from (d) of another 10 % at low and high p 2 . However, in the integrated width the latter 
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averages out almost completely. The uncertainties from (a) to (d) on / + + /~ are of 
comparable size, while the effect of radiative corrections is still unknown. 

The sum rules fl97[) and (|98|) for the B — * n form factors are formally converted into sum 
rules for the D — > tt form factors by replacing b with c and B with D. The input parameters 
are taken over from the calculation of fp in section 2. In addition, one has to rescale the 
wave functions from pb ~ 2.4 GeV to p c ~ 1.3 GeV as specified in Appendix 1. The allowed 
Borel mass window is 3 GeV 2 < M 2 < 5 GeV 2 . With this choice, the D — > n form factor at 



p = is predicted to be [42], [43 



/+(0) = 0.68. 



(101) 



0(a s ) corrections are not included here. The momentum dependence of / + in the range 
< p 2 < m 2 — 0(1 GeV 2 ) is shown in the next section, together with an extrapolation to 
higher p 2 . 

Other important applications of light-cone sum rules include the estimate |4l|] of the 
B — > K form factor which determines the factorizable part of the B — > J /%pK amplitude, the 
prediction [34] of the matrix element of the electromagnetic penguin operator for B — > K*^, 
and the more recent calculation ]33] of the B — > p form factors. In the latter work, essentially 



the same procedure is applied as outlined above. However, the relevant weak currents and 
light-cone wave functions are very different. For illustration, one has to deal with the 



following matrix elements |33| , [34], |37| 



(p + (q, A) | u(0)a^d(x) | 0) = -ift{e(\)pq v - e(A),g M ) / du e l ^0 ± (n) 



(102) 



(p+(q,X)\u(0) ltl d(x)\0) 



qx q 

+ U\), - qj^) fpmpf du e^ x g { l\u) , (103) 



(p + ( 9 ,A)|«(0) 7 ^d(z)|0) 



e(\yq p x°f p m p I du e^ x g [ l\u) 



(104) 



Here, 0j_ and 0n are twist 2 wave functions of transversely and longitudinally polarized 

p mesons, respectively, while gf 1 and g^ are associated with both twist 2 and twist 3 
operators. Higher-twist components and three-particle wave functions are still missing, as 
are perturbative corrections. 

There are four independent B — > p form factors: 



{p(X)\u^(l - l5 )b\B) 



-i(m B + m p )Ax{t)e{\)^ + 



iA 2 (t) 



(e(A)p B )(ps + Pp) 



iA 3 (t) 



Tfi B + m„ t/ " 



{z{>)Pb){Pb-Pp)l 



(105) 



2V(t) 



£{X)aPBpP P1 
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with t = (pb — Pp) 2 - Only V(t), A 1 (t) and A 2 {t) contribute to the semileptonic decay 
B — > plvi for I = e,fi. These form factors are shown in Fig. 11. They all are predicted 
to rise with momentum transfer t, contrary to some earlier claims based on three-point 
sum rules. This contradiction was studied in detail in |33j with the conclusion that the 
three-point sum rules are not reliable for reasons indicated in section 4. 



6 B*Bn and D*Dtt couplings 

The strong couplings of B and D mesons to pions have been studied with different variants 
of sum rules and in a variety of quark models. We shall use the following definition of the 
B*Bn coupling constant: 

(B*°n-\B-) = -g B . Bv (q.e). (106) 
The couplings of the different charge states are related by isospin symmetry: 

9B* Bn = 9b*°B-tt+ — —V^9B*°B°Tr = ^9B*-B—k° = —9B*-B°tv- ■ (107) 



In [ji2| , a light-cone sum rule for gB*Bn has been suggested which is derived from the 
same correlation function ( |55| ) as the corresponding B — > tt form factors and which depends 
on the same nonperturbative input. The key idea is to write a double dispersion integral 
for the invariant function F(p 2 , (p + q) 2 ). Inserting in ( |5BD complete sets of intermediate 
hadronic states carrying B and B* quantum numbers, respectively, and using the matrix 
elements @, (f41~D and ( |106|) , one obtains 

F( 2 ( + ) 2 )= m \ m B*fBfB*9B*B^ r p h (s 1 ,s 2 )ds 1 ds 2 , 10g , 

m b (p 2 - m 2 B ,)({p + q) 2 - m B ) 7s (s x - p 2 )(s 2 - (p + q) 2 ) 

+ {subtractions) . 

Again, the first term is the ground-state contribution and contains the B*Btt coupling, 
while the spectral function p h (si, s 2 ) represents higher resonances and continuum states 
in the B* and B channels. The integration region in the (s\,s 2 ) - plane is denoted by 
E. The subtraction terms are polynomials in p 2 and/or (p + q) 2 which vanish by Borel 
transformation of ( |108j ) with respect to both variables p 2 and (p + q) 2 . The transformed 
hadronic representation of F is given by 



F(M 2 ,M|) = B M? B Mi F(p 2 ,(p + q) 2 ) 



m 2 B m B *fBfB*9B*Bn 



cxp 



m 



B* 



M 2 



TTTb 

M 2 



+ 



cxp 



Si 



Ml 



S2 



P h (si, s 2 )ds±ds 2 



(109) 



where M\ and M| are the Borel parameters associated with p 2 and (p + q) 2 , respectively. 
The same transformation is also applied to Q9l)| ) yielding Fqcd{M 2 , M|). Then, the sum 
rule for 9b*Bw results from the equality of F and Fq C d, and continuum subtraction. This 
last step deserves some more explanations. 

Let us consider the twist 2 expression for F given in (|64]). In order to write it in the 
form of a double dispersion relation, 

F Q cn(p 2 , (P + qf) = f~ r - t\ ^ P QCD ^ s 2 ) , (110) 



ml Si -p 2 Jm\ S 2 - (p + q) 2 
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Figure 11: The B —>■ p form factors: predictions from light-cone sum rules (solid) in com- 
parison to lattice results j73|/ . The dashed curves indicate the uncertainties in the sum rule 
results (from ffijjj). 
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we change u to s = (ml — p 2 )/u + p 2 , and get 



Fqcd(p 2 , (p + qf) = rribf^ 



ds <^7r(s) 



Iml (s~p 2 )(s - (p + q) 2 ) ' 
In general, the wave function (p^{u) can be expressed as a power series in (1 — u): 



iir 



u 



s — mi 



s — p z 



(112) 



Substituting ( |112j) in (|111|) and introducing formally two variables s\ and s 2 instead of s, 
one reproduces the double integral ( |110[ ) with 



lfa k 



Now we Borel-transform ( |110| ), 

F QCD (M 2 J M 2 ) = m b f n J2 d 8l 



( Sl -m 2 b ) k 6^( Sl -s 2 ) 



(113) 



r(jfc + i) 



x 5 (k \ Sl -s 2 ) 



cxp 



_£l «2_ 

Ml Ml 



(114) 



introduce again new variables s = s\ + s 2 and v = Si/s, and use the <5-function to integrate 
over v. The result is 



F QCD (Ml Ml) = m b U E J 2 J S 



dv 



mi 



v — 



x exp 



svM 2 + s(l - v)M 2 
M'-Mq 



(115) 



v=l/2 



For Ml = Mf = 2M 2 the u-dependence of the exponent disappears, and the differentiation 
of the bracket gives a factor k\. As a consequence, ( |115|) reduces to 



F QCD (M 2 J M 2 )=m b f 7T J2 



2 k+i 



2mi 



dse = m b f 7T Lp 7T (u )M 2 e ^ 



;ii6) 



with Mo = 1/2. For arbitrary values of Ml and Mf a similar expression is obtained, with 
u = Ml I (Ml + Mf) and M 2 = Mf Mf /(Mf + Mf). 

With these manipulations, and the replacement of the integral over p h in ( |109|) by a 
corresponding integral over p® CD , it is straightforward to subtract the contributions from 
excited and continuum states from ( [L14f) . The remaining integral is restricted to the region 
below a given boundary in (si, s 2 ). For the latter one may take 



(si) a + (s 2 y < ( So y. 



(117) 



For a = 1, the duality region is a triangle in the (si, s 2 ) - plane, while for a — > oo it is a 
square. Since the spectral density ( |113j ) vanishes everywhere except at s± = s 2 , it is actually 
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irrelevant which form of the boundary one chooses provided the length of the duality interval 
at S\ = S2 is the same. For example, the triangle with so = is equivalent to the square 
with sq = Sq . Here, we take s£ as the effective threshold in both the B and B* channels. 



Repeating the steps following ( p. 14 ) one obtains an expression similar to (|115 ), but with the 
upper limit of integration in s lowered to 2s$ and with the addition of surface terms. The 
latter disappear for Mf = M 2 , in which case the subtracted invariant function replacing 
( |116|) is simply given by 



F subtr (M 2 ,M 2 ) 



2 fc+i 



2s B 



2m\ 



ds e 



m b f w ip^(u )M 2 



-4 



:ii8i 



However, in general, the proportionality of Fqcd{M\, Mf) to the wave function <p n at u = 
Mf/(Mf + Mf) found in (|116|) is destroyed by continuum subtraction. For this reason we 



adopt the particular choice Mf 



Mf 



[uq = 1/2) in what follows. 



So far we have solved the problem only for the twist 2 component of the sum rule. 
Unfortunately, the subtraction procedure explained above does not work for terms in Fqcd 
which contain higher powers of l/{m 2 — (p + uq) 2 ) (see fl9T|)). The main problem is that 
the corresponding double spectral densities are not concentrated near s\ = s 2 , making the 
continuum subtraction rather complicated in these cases. For further discussion we refer 
the reader to |50| . On the other hand, this difficulty only concerns higher- twist terms 
which contribute relatively little to the sum rule. Hence, to a good approximation one may 
disregard the continuum subtraction in these terms. 

Applying these recipes to Fqcd as given in ( p5| ) one ends up with the following light-cone 
sum rule for the B*Bn coupling: 



f bJ b*Qb*b IT 



m 



m B mB* 



-t9t 
— e m2 



+ e 
where 

and 



nib V <j b du I nibj n 



2 /3tt t g 



tf(«o) 



fx n m b 



ip a [uo) + 2u g 2 (u ) 



Ami 



no 



g 2 {v)dv ) +/f (u ) 



110=1/2 



da.\ 



(P3n(ai, 1 - u ,u - ax) f 1 " 011 , f^{ax, 1 - ax - a 3: a 3 ) 

CLCt 3 



/f(«o) 



Uq — OiX 

da 3 

dax / 

o JuQ—ai a 3 



«3 



[2<p ± (ati) - (f\\{ai) + 2(p 1 _{a i ) - <f\\{ai)] 



;ii9) 



(120) 



(121) 



are the contributions from the three-particle amplitude (|9"3"|). As a further simplification, 
G-parity implies g 2 (uo) = ( u o) — at u = 1/2, whence these terms can be dropped in 
the above sum rule. 



30 



The wave function f w (u ) at the symmetry point u = 1/2 also enters the sum rules for 
other hadronic couplings involving the pion. We take the value 



^(«=2) = 1.2±0.2 (122) 



obtained from the light-cone sum rule for the pion-nucleon coupling [3S]. This value is 



consistent with the choice of the coefficients 02,4 given in (AG). For the remaining parameters 
we use the same values as in the calculation of the form factor / + in section 5. With this 
choice the sum rule (|119|) yields 



f bJ b* Qb* b-k = 0.64 ± 0.06 GeV 2 , (123) 



9B* Bit = 29 ± 3 , (124) 



or 



where fs and /#* have been substituted from (p9|) and (|43|) , respectively. The uncertainties 
quoted above are to be interpreted as the range of values corresponding to the allowed 
window in the Borel mass, 6 < M 2 < 12 GeV 2 , in which the excited and continuum states 
contribute less than 30 % and the twist 4 corrections do not exceed 10 %. Also these 
predictions can and should be improved by calculating the radiative gluon corrections. 



The sum rule (|119|) for 9b*b-k is translated into a sum rule for 9d*d-k = 9d*+d°-k- by 
formally changing b to c, B to D, and B* to D* . With ( |122| ) and the input parameters from 
the calculation of the D —>■ ir form factor, one finds 

f D fD*9D*Dn = 0.51 ± 0.05 GeV 2 , (125) 

and using the LO estimates for f D and fo* from ( [4"0"D and fl4"4]), respectively, 

9d*Dtt = 12.5 i 1.0 . (126) 

The variation of the numerical results with M 2 in the allowed interval 2 < M 2 < 4 GeV 2 
is again quoted as an uncertainty. 

From fljgp one can calculate the width for the decay D* — > D7r. The predicted value, 

T(D* + ^ D°tt + ) = 9d * d I I q I 3 = 32±5keV, (127) 

lies well below the current experimental upper limit, 

T(D* + -> L)°7r + ) < 89 keV , (128) 

which is derived from the upper limit T tot (D* + ) < 131 keV and from the branching ratio 
BR{D* + — > D°7i + ) = (68.3 ± 1.4) % |5T]| . Predictions for other charge combinations are 
readily obtained from ( |127| ) and isospin relations analogous to ( p.07| ). Accounting also for 
the differences in phase space, one expects 

T(D* + -> D°n + ) = 2.2T(D* + -> D + tt°) = 1.44 r(D*° -> DV) . (129) 
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Table 2: Theoretical estimates of the strong B*Bir and D*Dir couplings (from H^J). 



Reference 
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9B*Btt 
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a QCD sum rules in external axial field or soft pion limit. 

* including perturbative correction to the heavy meson decay constants. 

b Quark model + chiral HQET. 

c Chiral HQET with experimental constraints on D* decays. 

d Extended NJL model + chiral HQET . 

e Quark Model + scaling relation. 

' Relativistic quark model. 

9 Bag model. 

" SU(4) symmetry. 

1 Reggeon quark-gluon string model. 

k Experimental limits 



32 




p 2 [GeV z ] 



Figure 12: The B — ► 7r /orm factor f + : direct sum rule prediction (solid), and single-pole 
approximation normalized by the sum rule estimate for gs*Bn (dashed). 



In Table 2, we summarize the numerical results discussed above and compare them with 
other estimates. One observes significant differences. Some of the predictions are rather 
close to the experimental upper limit, some even violate it. It would be very interesting to 
have more precise data. 

In contrast to gD*Dn, the coupling constant gB*B-w cannot be measured directly, since 
the corresponding decay B* — > Bit is kinematically forbidden. However, the B*Bir on-shell 
vertex is of great importance for the understanding of the behaviour of the B — ^ n form 
factors at large momentum transfer. Near the kinematic limit the form factor / + is expected 
to be dominated by the B* pole. The single-pole approximation given by 



f + ( P 2 



fB*9B*Bir 



2rriB*{l —p 2 /m 2 B ,) 



(130) 



is illustrated in Fig. 12 taking gB*Bn from (|124|) and /#* from fl43T). Extrapolation of the 
single-pole model to smaller p 2 matches quite well with the direct estimate from the light- 
cone sum rule fl97f ) at intermediate momentum transfer p 2 = 15 to 20 GeV 2 . This provides 
us with a consistent and complete theoretical description of / + . The extrapolation for the 
D — > 7i form factor using the analogous single-pole formula with go'D-w from ( |126|) and fo* 
from is shown in Fig. 13. Also in this case we find the direct sum rule result and the 
pole model to match nicely at p 2 ~ 0.7 GeV 2 . 

Finally, we refer the reader to for a similar calculation of scalar and axial B meson 
couplings yielding 



T(B(Q- 



B 



TV 



r(B(r 



B*n) ~ 360 MeV 



(131) 



and to p3| for an investigation of the B*Bp coupling. 



33 




p 2 [GeV z ] 



Figure 13: The D — > 7r /orm factor f + : direct sum rule prediction (solid), and single-pole 
approximation normalized by the sum rule estimate for Qd*d-k (dashed). 



7 Exclusive semileptonic decays, V u b, and all that 

With the B — > tt, D — > tt, and also B — + p form factors at hand, we are now in the position 
to predict the widths and differential distributions for exclusive semileptonic decays. In the 
case of f + (p 2 ), the sum rules (|97|) and ( |119| ) together with the single-pole approximation 
( |130| ) provide a complete description as illustrated in Fig. 12 and 13. For convenience, we 
have fitted parametrizations of the form 



f + (p 2 



/ + (0) 



ap 2 /m 2 P + bp^/rrip 



(132) 



mp denoting the B and D meson mass, to the theoretical results plotted in these figures. 
For the B — > 7r form factor we get 



/ + (0) = 0.27, a = 1.50, b = 0.52. 



(133) 



Here, the NLO correction to the twist 2 contribution has been included, and / + (0) has 
been fixed at the value given in (|100|) . At p 2 < 17 GeV 2 the fit reproduces the prediction 
of the light-cone sum rule (|97f) , while at p 2 > 20 GeV 2 it coincides with the single-pole 
approximation ( |130j ). In Fig. 14 the interpolation ( |132| ) is shown in comparison with recent 
lattice results [J The agreement is very encouraging. The analogous fit of (|132|) to the LO 
D — ► 7r form factor plotted in Fig. 13 yields 



/+(0) = 0.68, a = 1.16, 6 = 0.32 



(134) 



'For a comprehensive review see the article by J.M. Flynn and C.T. Sachraida in this volume, ref. 
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where / + (0) has been kept fixed at the value given in ( |101|) . 
The distribution of the momentum transfer squared in B 



ttIvi is given by 



dr 

dp 2 



G 2 \V r 



ub 



{p 2 - mf) 2 ^E 2 - m 



24vr 3 



p A m 2 B 




1 + m^Ei-mi) f + (p 



+ 



3m? 



m 



2\2 



ip A 




(135) 



with E n = (m 2 B + m 2 — p 2 )/2rriB being the pion energy in the B rest frame, and the 
form factors being as defined in (^) and (|90D. For / = e or fi, the form factor /° plays a 
negligible role because of the smallness of the electron and muon masses. Another interesting 
observable is the distribution of the charged lepton energy E\ in the B rest frame: 



dEi 



G 2 IK 



ub 
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.2 
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dp 2 
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4m B (p 2 + AEf) + 
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m B 



2m\ + p 2 



mi 
m B 

2m 2 



8m B Ei- 3p 2 + Ami 



ml 



4m B E[ + m 2 



m B 
f + (p 2 )f-(p 2 



f + ( P 2 



+ 



mt 



m B 



(P 2 



m. 



f-(p 2 



(136) 



with the integration limits pLax = m B (Ei ± J Ef — mf) + 0(m 2 ). The terms proportional 

min " 

to the pion mass squared not shown explicitly are taken into account in the numerical 
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calculations. Here, terms involving the form factor / are suppressed by the lepton mass if 
I — e, [A. 
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Figure 15: Distribution of the momentum transfer squared in B — > nlui (I = e,fi). 
dashed curves indicate the theoretical uncertainty discussed in detail in the text. 



The 



The above decay distributions are shown in Fig. 15 and 16 for I = e, /i. For the integrated 
width one obtains 

T(B -> nlut) = (7.5 ± 2.5) \V ub \ 2 ps" 1 , (137) 

where the theoretical uncertainty reflects the uncertainties in / + (added in quadrature) 
which have been discussed in detail in section 5.. 

Contrary to the semileptonic decays into e and //, the decay B — > nf u T is quite sensitive 
to the scalar form factor f°. Unfortunately, the single-pole approximation used to extrap- 
olate the sum rule result on / + to maximum p 2 cannot be applied to f°. This is because 
the scalar B ground state which should be about 500 MeV heavier than the pseudoscalar 
B lies too far above the kinematical endpoint p 2 = (m^ — m^) 2 of the B — > tt transition in 
order to dominate the form factor. Nearby excited resonances and nonresonant states are 
expected to give comparable contributions. For illustrative purposes J|3|], we extrapolate 
the form factor f° linearly from the maximal value p 2 = 15 GeV 2 at which the sum rules 
) still hold to the value at p 2 ~ m 2 B dictated by the Callan-Treiman limit [69| : 



and 



lim n f{p 2 ) = y- = 1.1 to 1.6 



'138) 



where we have used the conservative estimate fs = 150 to 210 MeV. This rough extrapo- 
lation of the sum rule result is plotted in Fig. 17. Also shown are lattice data. They are 
systematically lower than our expectation. It should however be noted that in contrast to 
/ + , NLO effects are still missing in the sum rule for f°. 



36 



Fig. 16 shows the resulting distribution of E T . As anticipated, the spectrum is quite 
sensitive to the large p 2 behaviour of f°, and may thus allow to determine, or at least 
constrain the scalar form factor experimentally. The integrated partial width is 

F(B -> txtv t ) = (6.1 ± 0.4) \V ub \ 2 ps~\ (139) 



yielding, together with (|137| ), the ratio 

T(B — ► 7TTZ/ T 



T(B — > 7rev e ) 



0.75 to 0.85 . (140) 



This ratio is independent of V u b, and less sensitive to uncertainties in the sum rule parameters 
than the widths themselves. The range quoted above corresponds to the variation of /° 
within the two extrapolations considered in Fig. 17. 




Figure 16: Distribution of the lepton energy in B — > -nlv\ for I = e,/i (solid) and I = t 
(dashed). In the latter case the two curves correspond to the two extrapolations of f° shown 
in Fig. 17. 



Recently, the CLEO collaboration |7(J has reported the first observation of the semilep- 
tonic decays B — > irlv and B — > plv (I = e, /x). From the measured branching fraction 
BR(B° -> Tt'l+Pi) = (1.8 ± 0.4 ± 0.3 ± 0.2) ■ 10" 4 , and the world average of the B° lifetime 
[pTf| , r B o = 1.56 ± 0.06 ps, one derives 

T(B° -> n-l + ui) = (1.15 ± 0.35) • 10" 4 ps" 1 , (141) 

where the errors have been added in quadrature. Comparison of (|141|) with (|137|) yields 

|y u6 | = 0.0039 ± 0.0006 ± 0.0006. (142) 
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Figure 17: TTie B ^ n form factor f°: direct sum rule estimate (solid) and linear extrapo- 
lations to the limit ( \138j ) (dashed). The lattice results are from jI3|/ . 



Here, the first (second) error corresponds to the current experimental (theoretical) uncer- 
tainty. 

A similar analysis can be performed for B — > plv\. Fig. 18 shows the distribution of the 
momentum transfer squared obtained from the form factors plotted in Fig. 11 |33]. The 
integrated width predicted in is 

T(B -> plv { ) = (13.5 ± 4) \V ub \ 2 ps- 1 . (143) 

Comparison with the CLEO result, 

T(B° -> p~l + vi) = (1.60 ± 0.6) ■ 10~ 4 ps" 1 , (144) 

derived from the measured branching ratio BR(B° — > p~l + i>i) = (2.5 ± QA~t Q ' 7 ± 0.5) • 
10 -4 and the B° lifetime also used in ( |141| ), gives 

\V ub \ = 0.0034 ± 0.0006 ± 0.0005. (145) 

Within errors, the values of \V u b\ extracted from the two exclusive semileptonic decays are 
nicely consistent with each other, and also coincide with the inclusive determination of V u b- 
A further test of the consistency between theory and experiment is provided by compar- 
ing the ratio (/ = e, v) 

T(B° -> p-l+vi) , , 

calculated from ( |143| ) and ( |137| ) with the ratio 

BR(B° - p~l+u e ) 1 



BR{B° -> n-l+u e 



1.4+^ ±0.3 ±0.4 (147) 
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observed by CLEO. Although the test is passed by the present estimate, the uncertainties 
on both sides are still too big to draw any firm conclusion. 
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Figure 18: Distribution of the momentum transfer squared in B 
show the theoretical uncertainty (from fiSBjJ). 



pev e . The dashed curves 



Turning to semileptonic D decays we show, in Fig. 19, the distribution of the momentum 
transfer squared in D — > irlui. For the integrated width we get 



T(D -> nlui) = 0.16 \V cd \ 2 ps" 1 = 8.0 • 10~ 3 ps' 



:i48i 



using \ V c d\ = 0.224 ±0.016 [51]. This prediction should be compared with the experimental 
result 



T(D° 



7T e u P 



(9.2t|l) • 10- 3 ps- 



derived from the branching ratio BR(D^ 



(149) 

TT-e+u) = (3.8j£g) • 10~ 3 and the lifetime 
r D o = 0.415 ± 0.004 ps |5T|. The present theoretical uncertainty is estimated to be of the 
order of the experimental error. In other words, the CKM-suppressed exclusive semileptonic 
D decays are not yet measured precisely enough to really challenge theory. In addition, there 
is a further demand for better data. One may use a precise measurement of D — ► nlui to 
constrain the light-cone sum rule for the D — * 7r form factor. Since the evolution of the pion 
wave function and other scale-dependent input quantities from the charm to the bottom 
scale is well under control, this is a promising way to considerably improve the sum rule 
predictions on B — > 7r transitions. 



8 Heavy quark limit 

The light-cone sum rules described in section 5 provide a unique possibility to investigate the 
heavy-mass dependence of the B — > tt form factors. To this end, one employs the following 
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Figure 19: Distributions of the momentum transfer squared in D — > irlvi (I = e,fi). 



scaling relations for mass parameters and decay constants: 

tub = m b + A , Sq = ml + 2m b ujQ , M 2 = 2m b T , (150) 
fB = fB/y/rnb, /s* = fB*/\frn b , (151) 



where in the heavy quark limit A, u , t, f B and /#* are m^-independent quantities. With 
these substitutions it is straightforward to expand the sum rules (|97|) and ( j9~8D in m&. In both 
cases, the light-cone expansion in terms of wave functions with increasing twist is consistent 
with the heavy mass expansion, that is the higher-twist contributions either scale with the 
same power of mj as the leading-twist term, or they are suppressed by extra powers of m&. 

We find that the asymptotic scaling of the form factors differs sharply at small and 
large momentum transfer | 43| . At p 2 = 

/ + (0) = /°(0)~m 6 - 3/2 , (152) 

/+(0)+r(0)~m 6 - 3/2 , (153) 
whereas at p 2 = ml — 2m b x, X being independent of m b , 

tip 2 ) -ml 12 , (154) 
f+(p 2 ) + f-(p 2 ) ~ /V) ~ < V2 • (155) 

The sum rules thus nicely reproduce the asymptotic dependence of the form factors / + and 



/ on the heavy quark mass m b derived in [58, py| for small pion momentum in the rest 



5 see also ref. |Hfl 
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frame of the B meson. In addition, the sum rules allow to investigate the opposite region of 
maximum pion momentum where neither HQET nor the single-pole model can be trusted. 
In particular, from the sum rule point of view it is expected that the excited and continuum 
states become more and more important as p 2 — > 0. This is reflected in the change of the 
asymptotic mass dependence. Claims in the literature which differ from ( 152|) to ( |155| ) are 
often based on the pole model and therefore incorrect in our opinion. A similar analysis has 
been carried out in p3|| with essentially the same conclusions. 

Furthermore, the sum rule (|119|) for the B*Bn coupling constant suggests 



9B* Bit ~ m b 



(156) 



in agreement with the expectation from HQET [[57], |58, 7lfl . At the expense of additional 



parameters such as A, uq, etc. light-cone sum rules can also be used to determine the I/thq 
corrections. For gB*Bn and gD*Dn, a simple quantitative estimate of the latter is obtained 



by fitting the numerical results ( 124 ) and ( |126j ) to the form 



9B* Bit 



^■sfi + A) 

U V m B 



and the analogous expression for go'D-n- This yields |42] 

£ = 0.32 ±0.02, A = (0.7 ± 0.1) GeV 



(157) 



(158) 



In Table 2 , the above value of the reduced coupling constant g is compared to the results in 
other approaches. The next-to-leading term in the heavy mass is sizeable, increasing from 
about 15% for gB*Bw to 40% for gwo-K- 

Finally, in the heavy quark limit the ratio 



gB*BirfB*\/ m D 
gD*DnfD*\/ m B 



(159) 



is expected to approach unity. Moreover, it has been shown J72[ that r is subject to 1/rriQ 
corrections only in next-to-leading order. Indeed, the ratio r = 0.92 derived from the 
light-cone sum rules deviates surprisingly little from unity, in agreement with the HQET 
expectation. 



9 The nonfactorizable amplitude for B —> J/ipK 

As a final example for applications of QCD sum rules to exclusive heavy meson decays we 
consider nonleptonic two-body decays. This class of processes is theoretically much more 
complicated than the (semi)leptonic decays discussed so far. As compared to the latter, 
effects from (a) hard gluon exchange at short distances, (b) soft interactions of quarks and 
gluons including nonspectator effects, (c) hadronization, and (d) final state interactions 
among the hadronic decay products change things considerably. Up to now, only the hard- 
gluon effects can be systematically taken into account in the framework of improved QCD 
perturbation theory. The result is an effective weak Hamiltonian at the physical scale 
fi ~ itiq <C rriyy, given by a sum of local operators with renormalized Wilson coefficients []. 

6 See, e.g., ref. {§. 
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We restrict our discussion to the decay B — > J/ipK which will play an important role at 
future .B-factories and which brings the main theoretical difficulties to light. The piece of 
the effective Hamiltonian relevant for this decay mode may be written in the form 

H w = ^= Vcb V; s {{c2 + j)0 2 + 2 Cl 2 } , (160) 

where 

OM = (cT^c)(sT p b), <5 2 (/i) = (cT^c)(sT p ^b) (161) 

with T p = 7 P (1— 75). The Wilson coefficients Ci(p) contain the effects from QCD interactions 
at short distances below the scale set by the inverse 6-quark mass. The hadronic matrix 
elements of the four-quark operators ( |161| ) are supposed to incorporate the long-distance 
effects (b) to (d). The problem of calculating these matrix elements is extremely demanding 
and still far from a satisfactory solution. 

In a radical first approximation, one may factorize the matrix elements of Hy/ for B — ► 
J/ipK into products of hadronic matrix elements of the currents that compose Hw- Strong 
interactions at scales lower than /1 between quarks entering different currents as well as 
nonspectator effects are thereby completely neglected. Moreover, the matrix element of the 
operator O2 vanishes because of colour conservation so that 

(J/1>K I H w I B) = ^f cb V: s (c 2 {fi) + ^\ (J/4jK I OM I B) . (162) 

The factorized matrix element of the operator O2 is given by 

{J/ijjK I 2 (ji) I B) = (J/V I cT?c \ti){K\ sT p b \ B) = 2^/+^m^ ■ q) , (163) 
where 

U = 405 MeV (164) 

is the decay constant determined by the leptonic width F(J/ip— > l + l~) = 5.26 ± 0.37 keV, 
and 

f+^ K = 0.55 ± 0.05 (165) 
is the B — ► K form factor at the momentum transfer p 2 = estimated [O] from a light- 



cone sum rule similar to the one for the B — > tt form factor given in (|97|) . Obviously, 
denotes the J/ip polarization vector, and q the K four-momentum. 

Already at this point one encounters a principal problem: since the matrix elements of 
quark currents in ( |163| ) are scale-independent, the //-dependence of {J/ipK \ 02(/x) | B) 
which is supposed to cancel the //-dependence of the Wilson coefficients in ( |162j ) in order 
to give a physically sensible result, is lost. Hence, the above approximation can at best 
be valid at a particular value of fx which could be called the factorization scale fip. The 
conventional assumption 6 is fip = 0{m b ). 

Using the next-to-leading order coefficients Ci^ip) in the HV scheme with A^L = 225 



MeV from [73| and taking fi = m b — 5 GeV, one has 



c 2 (/i) + = 0.155 . (166) 
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Together with ( |164| ) and (|165|) , this yields 



BR(B -> J/ijK) = 0.025% , (167) 



a branching ratio which is considerably smaller than the measurements pi] , |74 

BR(B- -> J/V>#~) = (0.101 ± 0.014)% , (168) 

5^(5° -> J/^if ) = (0.075 ± 0.021)% . (169) 

The quantitative failure and the scale problem pointed out above imply that naive fac- 
torization of matrix elements does not work. Factorization has to be accompanied by a rein- 
terpretation of the Wilson coefficients. For decays such as B — > J/ipK, the short- distance 
coefficient C2(/x) + ci(/i)/3 is substituted by an effective coefficient a 2 which is supposed to 
incorporate possible nonfactorizable contributions. Phenomenologically p6| , a-i is treated 
as a free parameter to be determined from experiment. From ( |162| ), ( |163j ), and ( |168| ), the 
most precise of the two measurements, one finds 

| | =0.31 ±0.02 , (170) 

where the quoted error is purely experimental. The sign of a B remains undetermined. 
The above value is close to the outcome of a comprehensive analysis [|J] [] of nonleptonic 
two-body .B-decays. 

The big difference between the short-distance and effective coefficient ( |166| ) and ( |170[ ), 



respectively, points at the existence of sizeable nonfactorizable contributions. The latter are 



also needed to cancel or at least soften the strong /z-dependence of ( [166] ). A deeper study 
shows that the dominant nonfactorizable effects should arise from the matrix element of the 
operator 2 - Writing the latter in the convenient parametrization 

(J/il>K | 6 2 (p) | B) = 2fJ B ^ K { l i)m^ ■ q) , (171) 



and adding it to the factorized matrix element (|163|) of 2 , one gets for the complete decay 
amplitude: 

(J/^K \H W \B) = V2G V cb V; s a^ K Uf^ K m4e^ ■ q) (172) 
with the effective coefficient Iff 



a^ K = cM + + 2 Cl (/,)%^ • (173) 

/ B—+K 

In Fig. 20, we show the partial width for B — > J/ipK as a function of the parameter fstpx 
associated with (O2). Note that fsipK = corresponds to naive factorization, while the 
fitted value of a 2 given in (|170| ) implies 

IbMl 1 = m b) = +0-04 or - 0.12 . (174) 
For a smaller scale the value of Jb^k is slightly shifted to the right on the real axis, e.g., 

IbMv = \ m b) = +0-06 or - 0.09 . (175) 
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T(B -> J/ipK) [10-V 
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Figure 20: The partial width for B — > J/ipK as a function of fsipK parameterizing the 
nonfactorizable contribution to the decay amplitude. The horizontal dotted lines represent 
the average of the experimental widths given in the text. The solid (dashed) curves show 
the theoretical expectation for // = rrib/2 (fi = nib). The arrows indicate the QCD sum rule 
estimate. 



We see that a nonfactorizable amplitude of 10 to 20 % of the size of the factorizable one is 
sufficient to conciliate expectation with experiment. 

The theoretical calculation of a 2 and of the analogous coefficients for other two-body 
decays ]2B| is one of the most important tasks in heavy flavour physics. As a first step 



in this direction, we have undertaken a rough estimate of fs^K using again QCD sum 
rule methods 0. Following the general idea put forward in |f76| , we choose the four-point 
correlation function 

BW(P, q) = J d*x d 4 y d'z e ^-Mw<0 | T{j«(x)ji(y)6 2 (z)j*(0)} | 0) (176) 

with j,f = bi-f 5 u, jf = c-f u c, and & = u^^s being the generating currents of the mesons 
participating in the decay B — > J/ipK, and p+q, p and q being the respective four-momenta. 

Generalizing the procedure applied in section 2 and 3 to the two- and three-point cor- 
relation functions (^) and (pEBP, respectively, one writes a dispersion relation for ( |176| ) in 



terms of intermediate hadronic states in the B, J/ip, and K channel. The ground state 

7 We refer here to the so-called new model for the heavy-to-light form factors. 

8 This work was done in collaboration with B. Lampe. Results were already reported in |fT5|. 
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contribution to 

ft ( . (Q 1 3% I gXg I it I jMg I Q 2 1 gKg I jg I o) 

VlP ' gJ ' (m^-g2)( m 2 -p2)(m|-(p + 9 ) 2 ) + > l A "J 

contains the matrix element ( p. 71] ) of interest. In addition, one has contributions from excited 
resonances and continuum states denoted above by ellipses which lead to a very complicated 
singularity structure. On the other hand, in the Euclidean region Q 2 = —q 2 > 0(1 GeV 2 ), 
p 2 < 0, (p+q) 2 < one can expand the correlation function (|176|) in terms of local operators: 

r»,<?) = E^r(p.?^)(W> • (178) 

d 

Because of the restriction to short distances it suffices to keep only the operators with 
low dimensions given in fllBD . We have included all relevant operators up to d = 6. The 
corresponding coefficients C% u (p,q, fi) have been calculated from the diagrams shown in 
Fig. 21. By equating (|177|) and (|178|) one can derive a sum rule for the matrix element 
(J/iPK | 2 | B). 

There are two complications that are not present in the two- and three-point sum rules 
discussed above. One problem is the presence of a light continuum in the B channel below 
the pole of the ground state £?-meson. This contribution to (|177| ) can be associated with 
processes of the type B — > ll D*D s n — > J/ipK, where an intermediate state carrying D*D S 
quantum numbers rescatters into the final J/ipK state. Formally, in ( |176| ) it is created 
from the vacuum by the combined action of the operator product Oijf . As a reasonable 
solution we suggest to cancel this unwanted piece against those terms in the OPE ( |178 



with the quark content ccsq which develop a nonzero imaginary part at (p + q) 2 > 4m 2 . 
In the approximation considered this is the case for the four-quark condensate contribution 
represented by Fig. 21c. 

The second problem concerns the subtraction of contributions from excited resonances 
and continuum states in the remaining sum rule. A closer look at the J/i/j channel of the 
correlation function ( |176| ) reveals that the higher charmonium resonances contribute with 
alternating signs. Therefore, the usual subtraction procedure in which the dispersion integral 
over the excited and continuum states is approximated by its perturbative counterpart is not 
reliable here. In order to proceed we employ explicit, although rough models for the hadronic 
spectral functions. Since the number of additional parameters has to be manageable, only 
the first excited resonances are included in each channel besides the B, J/ip and K ground 



states. In total, in this approximation the correlator (|177|) contains three free parameters 
in addition to }'b^k- 

We then Borel transform (|177|) and (|178 ) in the 5-meson channel and take moments in 
the charmonium channel. In the i^-meson channel, q 2 is kept spacelike. Fitting the hadronic 



representation ( 177 ) to the OPE result (|178| ) for various values of the Borel mass M and 



q 2 , and for several moments, we find 

f mK = -(0.045 to 0.075). (179) 
The scale \x implicit in this estimate is set by the Borel mass. The central value in the 



allowed range turns out to be given by M ~ i/m^ — m\ ~ ~ 2.4 GeV. Substituting 



(|179|) in ( |173|) , and evaluating the short-distance coefficients c\^{\l) also at \i = M, one gets 
a B4>K = _ Q 2g + 3g _ ( ig tQ 31 ^ = _( _ 10 tQ o 22) , (180) 
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where the three terms in the first relation refer to the three terms in ( |173| ) in the same 
order. Interestingly, the sum rule approach seems to favour the negative solution for fs^K- 
Although in comparison with (|175|) our estimate falls somewhat short, the gap between 
theory and experiment is narrowed considerably as can be seen from Fig. 20. 





Figure 21: Diagrams determining the Wilson coefficients of the OPE of the correlation 
function IT^. The symbols are as in Fig. 1. 



Several comments are in order. Firstly, the nonfactorizable matrix element (171) is small 
as compared to the factorizable one given in ( |163| ), numerically, \fsipK I ' fB^K( m %)\ — 0-1- 
Nevertheless, it has a strong quantitative impact on because of the large coefficient 

|2ci/(c 2 + Ci/3)| ~ 20 to 30. Secondly, the factorizable amplitude proportional to Ci/3 
and the nonfactorizable one proportional to fs^K are opposite in sign and hence tend to 
cancel. In fact, if is taken at the upper end of the estimated range, the cancellation 

is almost complete resulting in a considerable enhancement of the branching ratio (|167|) . 
Note that both terms are nonleading in 1/N C . This is exactly the scenario anticipated 
by the 1/N C - rule for D decays JT7| which has found theoretical support by the global 
sum rule analysis of [7B|. In |78| , |T5f , a similar trend was claimed for B — ► Dti. Thirdly, 
our estimate yields a negative overall sign for a^ K in contradiction to naive factorization 
(first two terms in ( |173| )), and also to a global fit of the factorized decay amplitudes with 
two universal coefficients a\ and a 2 to the data It should be stressed, however, that 
in this fit the positive sign of a 2 actually results from the channels B~ — > D°tt~ , D°p~ , 
D*°ir~, and D*°p~, and is then assigned also to the J/ipK channel. This assignment may 
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not be correct. Certainly, the sum rule approach described above provides no justification 
for such an assumption. On the contrary, diagrams of the kind shown in Fig. 21 suggest 
some channel-dependence of the nonfactorizable matrix elements. For class II processes 
involving a<i = c<i + c\ /3 + 2c±f '/ f + (see ( |173j )), the channel-dependence is enhanced by the 
large coefficient 2c±, while the factorized matrix elements come with the small coefficient 
C2 + ci/3. A concrete numerical example is provided by ( |180| ). The opposite is the case 
for class I processes involving a\ = c\ + C2/3 + 2c2/// + : the nonfactorizable contributions 
are damped by the small coefficient C2, while the factorized term has a large coefficient. 
Therefore, a\ is indeed expected to be universal to a good approximation, but a 2 should 
exhibit some channel-dependence, in particular when comparing decays with very different 
final states such Dn and J/ipK. 

As a last remark, from the sum rule point of view there is also no simple relation 
between B and D decays. In H, arguments are presented which support a change of sign 
in ci2 when going from D to B decays. However, it is not obvious that these arguments 
hold independently of the particle composition of the final state. In the sum rule approach, 
for example, one has significant differences in the OPE of the correlation functions such as 
(|1 76|) for B — > Dtt, B — > J/ipK, and D — > Kit as can be imagined from Fig. 21. From this 
point of view at least, the relation of B and D decays is expected to involve more than just 
a change of mass scales. 



10 Conclusion 

QCD sum rule techniques have proved to be very useful in calculating hadronic matrix 
elements for exclusive decays of B and D mesons. In this review, applications are discussed 
to decay constants, form factors, and amplitudes of nonleptonic two-body decays, as well 
as to the B*Bir and D*Dir couplings. We have not considered exclusive radiative decays 
such as B — > K*^j and B — > p7, but we at least want to mention that here sum rules have 
been employed to determine the matrix elements of the leading magnetic penguin operator 



[p0| , 131] ? |32| , |34[ and to estimate long-distance effects |80], fH], p2| . The above examples are 
by no means exhaustive. 

Using the sum rule results we have presented predictions on decay distributions and 
integrated widths for B — > nlui, B — > plui, and D — > ttIui. Comparison with the CLEO 



measurements [F0| of exclusive semileptonic B° decays yields values for V u b in good agree- 



ment with each other and with the determination from inclusive data H&jfll - Agreement 
between expectation and measurement is also found for the Cabibbo-suppressed semilep- 
tonic D° decay |5T| . 



From the sum rule estimate of the D*Dn coupling we have calculated the decay width 
for D* — > Dn and compared our result with other estimates. The experimental upper limit 
[pTf| is still about three times larger than the expected width. 

Furthermore, we have described an attempt to estimate nonfactorizable matrix elements 
of nonleptonic decay amplitudes with the help of sum rule techniques. Using B — > J/i/jK as 
a prototype example, it is argued that nonfactorizable effects play an essential role. With 
the sum rule estimate of the nonfactorizable contributions included the effective coefficient 
a^ K is found to be consistent with the value extracted from the experimental branching 
ratio [[STJ. However, the sign of predicted by the sum rule analysis is opposite to the 
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sign determined from data, if channel- independent, universal coefficients are assumed 
[0. We have explained why in the sum rule approach universality can be expected for a 1; 
but not for a<i- This issue certainly requires further clarification. As long as one is not 
able to actually calculate the effective coefficients, one cannot claim a complete theoretical 
understanding of the exclusive nonleptonic decays. 

Among the different variants of sum rules we have put particular emphasis on the light- 
cone sum rules which provide a powerful tool in problems involving a pion, kaon or p meson 
such as heavy-to-light form factors and couplings. In this approach, the light hadrons are 
described on mass-shell by a set of wave functions with different twist and quark-gluon 
multiplicity, representing distribution amplitudes in the fraction of the hadron light-cone 
momentum carried by the constituents. This avoids the notorious model-dependence of 
extrapolations from Euclidean to physical momenta in light channels. The light-cone sum 
rules also seem to be fully consistent with the heavy mass expansion in contrast to some 
of the sum rules based on short- distance expansion. Moreover, the derivation of light-cone 
sum rules is often technically easier than conventional sum rule calculations. 

We have addressed in some detail the present theoretical uncertainties, and the prospects 
for improvement. On the theoretical side, we see room for it by determining the nonasymp- 
totic features of the light-cone wave functions more accurately, and by including higher-order 
perturbative effects in the sum rules. Work is under way in both directions: the re-analysis 



of the twist 2 pion wave function in [S5|, the determination of p wave functions in |pq| , and 



the calculations of the 0(a s ) correction to the B — > ir form factor / + in |47|, f|8| are recent 
examples. On the experimental side, the advent of new and more precise measurements at 
future B and tau-charm factories should allow to tightly constrain the input parameters and 
to test the reliability of the sum rule approach in a very challenging way. Because of the uni- 
versality of the nonperturbative input it appears conceivable to decrease the uncertainties 
from presently 20 to 30 % to about 10 %. 

Finally, we have pointed out the encouraging agreement of lattice and sum rule calcula- 
tions in the case of fs, Jd, an d the B — > tt and B — > p form factors. This agreement should 
be enough motivation to join efforts. We believe that it would be very fruitful to combine 
the flexibility of the sum rule method with the rigorous nature of the lattice approach. 
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Appendix 1 

Here, we collect the formulae for the light-cone wave functions of the pion and specify the 
parameters. It is important to note that the asymptotic form of these functions and the 
scale dependence are given by perturbative QCD [[J7|, |4lJ . 

The twist 2 wave function <p n is expressed as an expansion in Gegenbauer polynomials: 



<p n (u, p) = Qu(l - u) 1 + a 2 {p)Cl /2 (2u - 1) + a A {p)Ct'\2u - 1) + 



,3/2, 



(Al) 
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where 



C 3 4 /2 (2u - 1) = — [21(2u - l) 4 - U(2u - If + 1] . 



C 2 3/2 (2 M -l) = |[5(2n-l) 2 -l] 
15 



The normalization is such that 



du(p w (u, /i) = 1 



(A2) 



(A3) 



The nonperturbative effects are contained in the coefficients a n . In LO, they are multiplica- 
tively renormalizable and have the following scale dependence: 



a n (n) = a n (n ) 



ln/b 



(A4) 



where 6=11 — 2rif/3 is the LO coefficient of the QCD beta function, nj being the number 
of active flavours, and 



In = Cf 



-3 



(n + l)(n + 2) 




(A5) 



are the anomalous dimensions |37|| . We see that a n (/i), n > 2 vanishes for fi — > oo. Therefore, 
these terms describe nonasymptotic features of the wave function ( |A1|) . 

The initial values of the nonasymptotic coefficients can be estimated from two-point 
sum rules [|37| for the moments / u n (p^{u^ jj)du at low n. The nonperturbative information 
encoded in the quark and gluon condensates is thereby transmuted into the long-distance 
properties of the wave function. Alternatively, one can determine the coefficients directly 
from light-cone sum rules for known hadronic quantities such as the irNN and Lupir couplings. 
For the numerical results shown in this review we have used the following estimate at 
Ho = 0.5 GeV @: 

2 

c^Oo) = 77, 04(^0) = 0.43 . (A6) 



On the basis of the approximate conformal symmetry of QCD it has been shown ||6| that 
the expansion ( |AT| ) converges sufficiently fast so that the terms with n > 4 are negligible. 

The scale \x to be used in the light-cone sum rules ( |97D and fl98|) for the B — > tt form fac- 
tors, and in the sum rule (|119| ) for the B*Btt coupling is somewhat ambiguous, in particular 
in LO approximation. As a reasonable choice, we take 



m™ — mi ~ 2.4 GeV . 



The analogous choice for the D meson form factors and coupling is 



fJL c 



m 



D 



mi 



1.3 GeV . 



(A7) 



(A? 



These scales characterize the typical virtuality of the b, respectively c quark, and coincide 
within a factor of two also with the respective value of the Borel mass M. With ( |A4| ) and 
l6T) , one gets for the evoluted coefficients 



02 (/•*!>) = 0.35, a 4 (/i 6 ) = 0.18 , 
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a 2 (/x c ) = 0.41, a 4 (/i c ) = 0.23 . (A9) 



The evolution of a„(/z) is also known in NLO, where mixing effects occur [|9|. Using the 
same input values (|A6|) and \i b from (|A7|), one finds 



a 2 (^ b ) = 0.22, a 4 (/i 6 ) = 0.084 . (A10) 

The twist 3 two-particle wave functions tp p and tp a are related to the three-particle wave 
function ip 3n by equation of motion. It is therefore sufficient to specify the latter |37, 45, 46 . 
Including the first three nonasymptotic terms, which is consistent with retaining 02,4 in </?„-, 
one has 

^(cti, a0 = 360aia 2 a 2 (l + cj li0 (/f)-(7ai3 - 3) 
+ ^2,o(A i )(2 - 4aia 2 - 8a 3 + 80:3) + LUi tl (fi)(3aia 2 - 2a 3 + 3a 2 ) (All) 

yielding 

<p p (u, n) = 1 + S 2 (m) o ( 3 (« " " X ) + ^4(m) o (35(« - m) 4 - 30(u - u) 2 + 3) (A12) 
2 8 

and 

<f a (u, 11) = 6uu\l + C 2 (fi)-(5(u - u) 2 - 1) + CM— (21(u - m) 4 - 14(u - w) 2 + 1)1 (A13) 
L 2 8 J 

with u = 1 — u and 

^2 = 30^^, S 4 = |-^-(4w 2 ,o-wi i i -2w 1)0 ), 

C 2 = A-(5 - ^ 1>0 ) , C 4 = l^(4c,2,o - • (M4) 
The parameter /W^) and the coefficients a** have again been estimated from sum rules 

m- 

/ 37r (lGeV) = 0.0035 GeV 2 , 
wi,o(lGeV) = -2.88, cj 2 ,o(lGeV) = 10.5, ^(lGeV) = . (A15) 
After renormalization p7| , f46|| to the relevant scales //& and /i c one has 

/ 37 r(At 6 ) = 0.0026 GeV 2 , u lfi (ji b ) = -2.18, oo 2 M = 8.12, = -2.59 , 

fz-K^c) = 0.0032 GeV 2 , u lfi (ji c ) = -2.63, cu 2i0 (/i c ) = 9.62, Wljl (/i c ) = -1.05 . (A16) 

The parameter /i^ = m 2 /(m u + m^) can be inferred from the PCAC relation (p2|): 

^(1 GeV) = 1.65 GeV, ^(^ c ) = 1.76 GeV, = 2-02 GeV . (A17) 

For convenience, we also list the complete set of twist 4 wave functions given in ref. (46j. 
It includes four three-particle wave functions specified by two parameters: 

(p±(ai,fi) = 305 2 {n)(a 1 -a 2 )al[^ + 2e{n)(l-2a 3 )] , 
<p\\(aii,fi) = l205 2 (/j)e(jj,)(ai - a 2 )aia 2 a 3 , 

£±K//) = 305 2 (/i)a 3 2 (l-a3)^ + 2£(/i)(l-2a3)] , 

<fi\\(ai,ii) = -1205 2 a 1 a 2 a 3 [^ + e{fi)(l-3a 3 )} , (A18) 
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and two two-particle wave functions related to the former by equations of motion: 



5 1 .6 

gi(u, fi) = -<5 2 (u)w 2 m 2 + -e(n)8 2 (fi)[uu(2 + 13mm) + 10m 3 lnw(2 — 3m + -u 

6 



+ 10u 3 lnu(2-3M + -M 2 )], 





10 



<5 2 (u)wm(m — m) . 



The parameter 5 2 is actually defined by the matrix element 

(7T\g s dG atl >y a u\0) = i5 2 f^ . 



Renormalizing the values |83j 



5 2 (lGeV) = 0.2 GeV 2 



and [39 



e(lGeV) = 0.5 

obtained from sum rule estimates to the relevant scales \i c and «b, one finds 

5 2 (u c ) = 0.19 GeV 2 , e(u c ) = 0.45 , 
5 2 (fi b ) = 0.17 GeV 2 , e(n b ) = 0.36 . 



(A19) 

(A20) 

(A21) 
(A22) 



(A23) 



Appendix 2 



Here, we give the explicit expressions for various subdominant contributions to the light- 
cone sum rule 



the surface term t A 



4m 2 



(m 2 — p 2 ) 2 



I +?LX) <7i(A) + 



Ami 



dgi(A) 



+ 



mi — p A 



M 2 

B ^2" 



i s o ~ P 2 ){jn 2 — p 2 ) du 



1 + ml +p* t + gp -p 
m\—p 2 \ M 2 



A ?(m 2 4- n 2 ) 1 

92(v)dv - 2 1 b J P B ] g 2 (A) , (Bl) 
i! [mi — p^){Sq — p z ) J 



the contribution of three-particle wave functions Jq 

f^p\M 2 ) = 



1 , r VcaQi^ + ua 3 - A) 

udu 1 , V2 

o J («i + ua^) z 



x exp 



m b — P 2 (l — «1 — M«3) 



(«! + ua 3 )M 2 



$ 3 (M,ai,M 2 ,p 2 ^ 



(B2) 
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where 



2/; 



3- 



fnm b 



1 - 



2 2 

mi - p 



(«! + ua 3 )M 2 



uM 2 



2<^±(ai) - ¥>||(o!i) + 2(p ± (oii) - {p\\{oti) 



(B3) 



and the perturbative correction 5 + (p 2 ,M 2 ) 

-, 1 "o , 2 

S + (p 2 ,M 2 ) = - / duip n (u,n) / — IrnTW— ,— ,w,p)exp — , (B4) 
7r J J mi K mi mi ' \M Z J 



where 



7T 



ImTi(ri,r 2 ,u,yu) = 5(1 - p) 



7T 2 -6 + 31n^ -2Li 2 (n) 
P 



+2Li 2 (l - r 2 ) -2 ( In ^— f^f + 2 (lnr 2 + i— ^ ) (2 ln(r 2 - 1) - ln(l - r x )) 



+o( P -i) 



1 — ri 
ln(p-l) 



p-1 



r 2 



1 _ _ . . . . . ml \ 1 



2 ( In r 2 + — - 2 - 2 ln(r 2 - 1) + In ^ . 
V r 2 p 2 / p-1 



r 2 — 1 



'(ri - r 2 )(p - n) 



lnp-21n(p-l) + l-ln 



m| 
P 2 



+2 



1 — ri 



(ri - r 2 )(r 2 - p) \ r 2 



ln^-21n^i^ -4^ + 2^- 







i r 
+ --- 




r 2 J 


p p 



+0(1 - p) 



2 lnr 2 + 21n(r 2 - 1) -ln-| 



m?\ 1 



r 2 



p 2 / p-1 



- 2 



1 — n 



(ri - r 2 )(r 2 - p) 



lnr 2 + 1 — 2 ln(r 2 — 1) — In 



P 2 , 



with 



p = ri + w (r 2 -r x ), f d Pf(p) Y~T^ = J 



1 l-r 2 
r 2 - p r 2 

1 



dp (/(p) -/(!)) 



1-p 



The surface term t ± appearing in the sum rule (j98[) is given by 



t ± (s*,p 2 ,M 2 )=ex V 



M 2 



(p a (A) 2m b g 2 (A) 



6 A p^m 2 — p 2 



(B5) 



(B6) 



(B7) 
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